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At  ft*  act 

A  computerized  Tidal  Current  Model  for  Long  Island  Sound  and  San  Francisco 
Fay  has  been  developed  by  the  Coast  Guard  Oceanographic  Unit.  The  model 
is  a  digitized  version  of  the  National  Ocean  Survey's  Tidal  Current  Tables 
and  the  Tidal  Current  Charts  for  Long  Island  Sound  and  Block  Island  Sound 
and  San  Orancisco  Bay.  Tests  were  conducted  to  verify  the  accuracy  of  these 
sources  of  tidal  current  information  using  a  drift  pole  survey  of  surface 
currents  for  Long  Island  Sound  published  by  the  Coast  and  Geodetic  Survey 
and  the  Army  Corps  of  Engineers  Hydraulic  Model  Bay  and  Delta  Model  for 
San  Francisco.  The  Long  Island  Sound  test  showed  complete  agreement  between 
the  model  and  observations  83.2%  of  the  time,  and  -100%  agreement  at  40% 
of  the  positions  tested.  The  San  Orancisco  Bay  test  showed  under  pre¬ 
diction  of  speeds  by  0.5  to  1.0  knots  in  3  of  4  regions  and  over  prediction 
by  -0.6  knots  in  the  4th  region.  The  predicted  directions  were  deflected  to 
the  right  of  the  observed  directions  75%  of  the  time  by  49.15  to  66.11 
degrees.  These  comaarisons  indicate  that  the  data  for  Long  Island  Sound 
is  reliable  and  that  for  San  Francisco  Bay  should  be  used  with  caution. 
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ABSTRACT 


A  computerized  Tidal  Current  Model  for  Long  Island  Sound  and  San  Francisco  Bay  has  been  developed  by 
the  Coast  Guard  Oceanographic  Unit.  The  model  is  a  digitized  version  of  the  National  Ocean  Survey’s 
Tidal  Current  Tables  and  the  Tidal  Current  Charts  for  Long  Island  Sound  and  Block  Island  Sound  and 
San  Francisco  Bay.  Test  were  conducted  to  verify  the  accuracy  of  these  sources  of  Tidal  Current  informa¬ 
tion  using  a  drift  pole  survey  of  surface  currents  for  Long  Island  Sound  published  by  the  Coast  and  Geodetic 
Survey  and  The  Army  Corp  of  Engineers  Hydraulic  Model  Bay  and  Delta  Model  for  San  Francisco.  The 
Long  Island  Sound  test  showed  complete  agreement  between  the  model  and  observations  83.2%  of  the 
time,  and  100%  agreement  at  40%  of  the  positions  tested.  The  San  Francisco  Bay  test  showed  under  predic¬ 
tion  of  speeds  by  0.5  to  1.0  knots  in  3  of  4  regions  and  over  prediction  by  0.6  knots  in  the  4th  region. 
The  predicted  directions  were  deflected  to  the  right  of  the  oberved  directions  75%  of  the  time  by  49.15  - 
66. 1 1°  .  These  comparisons  indicate  that  the  data  for  Long  Island  Sound  is  reliable  and  that  for  San 
Francisco  Bay  should  be  used  with  caution. 
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INTRODUCTION 


An  automated  system  for  providing  tidal  current  information  for  Long  Island  Sound  and  San  Francisco 
Bay  has  been  developed  for  use  in  the  Coast  Guard’s  Search  and  Rescue  (SAR)  planning  and  in  Marine 
Environmental  Protection  (MEP)  such  as  forecasting  pollutant  drift.  The  system  is  a  digitized  version  of 
the  National  Ocean  Survey’s  tidal  current  charts  for  Long  Island  Sound  and  San  Francisco  Bay.  The  new 
system  is  an  updated  version  of  the  original  system  developed  by  Morgan  et  al.  (1974)  for  Long  Island 
Sound  with  San  Francisco  Bay  added.  The  objectives  of  these  systems  is  to  provide  tidal  current  informa¬ 
tion  for  use  with  the  Computer  Assisted  Search  Planning  (CASP)  in  a  timely  manner.  The  current  values 
generated  are  the  tidal  currents  and  estuarine  flow.  They  do  not  include  currents  generated  by  any  other 
environmental  parameters. 

Long  Island  Sound  is  classified  by  Swanson  (1975)  as  having  a  semidiurnal  tide  with  a  slight  diurnal 
inequality.  The  inequality  in  the  range  is  reflected  in  the  magnitude  of  the  ebb  and  flood  tide  currents.  In 
San  Francisco  Bay  there  is  also  a  diurnal  inequality  which  is  more  pronounced  than  it  is  in  Long  Island 
Sound.  This  makes  for  a  greater  variation  in  the  magnitude  of  the  ebb  and  flood  tides  occurring  during  a 
single  tidal  day,  24  hours  and  50  minutes.  That  is  to  say  the  flood  tides  in  a  given  day  may  have  significantly 
different  magnitudes.  The  same,  of  course,  is  true  of  the  ebb  tides.  This  inequality  is  produced  primarily  by 
the  lunar  declination  and  is  independent  of  whether  the  moon  is  north  or  south  of  the  equator.  The  diurnal 
inequality  has  a  maximum  magnitude  at  approximately  14  days  intervals  (Defant  1961  and  Disney  et  al. 
1925).  This  diurnal  inequality  in  the  magnitude  of  tidal  currents  makes  long  term  drift  forecasting  a  very 
tedious  process.  It  is  necessary  to  update  drift  forecast  positions  at  intervals  of  one  hour  or  less  and  to  make 
multiple  runs  of  the  tidal  current  program  to  obtain  the  correct  tidal  vectors. 

The  average  period  for  the  flood  tide  in  Long  Island  Sound  is  5.9  hours  and  the  average  period  for  the  ebb 
is  6.6  hours  based  on  the  month  of  January  1975.  By  dividing  the  flood  tide  into  six  equal  periods  and  the 
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ebb  into  seven  equal  periods,  a  total  of  thirteen  tidal  hours  of  approximately  equal  length  can  be  produced. 
These  tidal  hours  are  0.98  and  0.94  solar  hours  in  length  for  the  flood  and  ebb  respectively.  As  a  result  the 
time  steps  may  be  treated  as  equivalent  to  solar  hours. 

The  flood  tide  in  San  Francisco  Bay  is  longer  than  the  ebb  tide,  6.71  hours  versus  5.71  hours  (Disney 
1925).  Even  with  this  inequality  of  ebb  to  flood  time  the  tidal  sequences  are  divided  into  equal  intervals  for 
both  sequences.  These  tidal  hours  are  roughly  equilavent  to  solar  hours  being  1.12  hours  for  flood  and  0.95 
hours  for  ebb. 

SYSTEM  DESCRIPTION 

The  tidal  currents  are  obtained  from  the  National  Ocean  Survey  tidal  current  charts:  San  Francisco  Bay 
6th  edition  1973  and  Long  Island  Sound  and  Block  Island  Sound  6th  edition  1977.  These  charts  are  based 
on  direct  current  measurements  with  Roberts  Current  Meters  and  Richardson  Current  Meters  (Hicks,  1967). 
The  horizontal  distribution  of  the  current  vectors  in  these  charts  is  too  random  to  permit  the  needed  pre¬ 
cision  required  in  Coast  Guard  SAR  and  MEP  applications.  To  overcome  this  weakness  current  values  were 
interpolated  between  known  values  and  a  more  complete  array  developed. 

The  resultant  current  vector  arrays  have  a  horizontal  resolution  of  2  x  2  nmi.  This  is  2’  of  latitude  by  3’ 
of  longitude  in  Long  Island  Sound  and  2’  of  latitude  by  2.5’  of  longitude  in  most  of  San  Francisco  Bay  and 
San  Pablo  Bay.  In  the  approaches  to  the  Golden  Gate  and  the  parts  of  the  Bay  near  Angel  Island,  Treasure 
Island.  Yerba  Buena  Island  and  Alcatraz  Island  the  horizontal  resolution  is  1  x  1  nmi  or  1’  of  latitude  by 
1.25’  of  longitude. 

To  use  the  models  an  operator  needs  certain  parameters  which  are  used  by  the  computer  model  to 
calculate  the  tidal  currents.  These  input  parameters  are  the  local  date,  time  and  position.  The  position  is 
given  to  the  nearest  minute  for  Long  Island  Sound  and  the  nearest  0.25  minute  for  San  Francisco  Bay.  The 
output  of  the  program  is  the  tidal  current  speed  and  direction  in  degrees  true  and  knots  plus  the  tidal  cur- 
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rents  at  the  reference  station  for  the  date  of  occurrence.  The  Long  Island  Sound  currents  are  referenced  to 
the  Race  and  the  San  Francisco  Bay  currents  are  referenced  to  the  Golden  Gate. 


The  computer  program  LITSF,  which  is  used  to  generate  the  tidal  currents,  is  a  computerized  version 
of  (a)  the  annual  Tidal  Current  Tables  published  by  the  National  Ocean  Survey  for  the  East  and  West  coast 
and  (b)  the  Tidal  Current  Charts,  also  published  by  the  National  Ocean  Survey.  The  main  part  of  the  pro¬ 
gram  is  the  National  Ocean  Survey’s  Tidal  Current  Prediction  Program  which  is  used  to  generate  the  tidal 
current  information  in  (a).  The  Tidal  Current  Charts  have  been  digitized  and  stored  on  a  disk  file  in  the 
computer’s  memory.  These  files  are  accessed  by  a  subroutine,  SLACK,  in  the  program. 

When  a  date,  time,  and  position  of  occurrence  is  fed  into  the  compute,  the  program  uses  the  date  to 
calculate  the  tidal  current  information  for  the  reference  station  in  the  Tidal  Current  Tables.  This  informa¬ 
tion  is  then  used  by  the  subroutine  SLACK  to  generate  the  correct  tidal  current  for  the  time  and  place  of 
occurrence.  In  short  the  computer  is  doing  the  same  steps  as  a  man  using  the  two  references  above  in  genera¬ 
ting  a  tidal  current  vector. 

A  listing  of  the  program,  digitized  tidal  current  charts  and  an  explanation  of  how  the  charts  are  digitized 
are  included  in  Appendix  III. 
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SPECIAL  FEATURES  OF  LONG  ISLAND  SOUND 

The  tides  in  Long  Island  Sound  are  co-oscillating  tides  displaying  a  mixture  of  standing  wave  and  pro¬ 
gressive  wave  characteristics.  The  loss  of  energy  due  to  bottom  friction  is  such  that  only  a  partial  reflection 
of  the  tidal  wave  occurs.  This  is  most  clearly  seen  in  Long  Island  Sound  in  the  increased  amplitude  of  the 
western  end  of  the  sound.  The  primary  tidal  current  and  tidal  waves  enter  the  eastern  end  of  Long  Island 
Sound  through  the  Race.  The  tide  belatedly  enters  the  Sound  from  the  west  by  way  of  New  York’s  Hell 
Gate.  This  tends  to  dampen  the  flood  and  ebb  currents  in  the  western  end  of  the  Sound.  As  the  tide  enters 
Long  Island  Sound  through  the  Race  frictional  forces  retard  its  propagation  through  the  area.  There  are 
very  strong  currents  in  the  area  of  the  Race  reaching  speeds  as  great  as  5.3  knots  (272.8  cm/s)  during  ebb  and 
4.3  knots  (221.3  cm/s)  during  flood.  High  water  requires  an  hour  to  travel  or  extend  from  Montauk  Point  to 
the  Race  and  three  and  a  half  hours  to  extend  into  the  relatively  small  Peconic  Bays  to  the  south  of  the 
Race.  The  lag  from  Montauk  Point  to  Throgs  Neck  is  three  hours.  As  a  result  the  flood  current  begins  at  the 
eastern  end  of  the  sound  and  moves  progressively  westward.  The  currents  at  any  given  point  in  the  Sound 
reach  their  maximum  midway  between  slack  water  times. 

The  high  current  speeds  in  the  region  of  the  Race  require  frequent  updating  of  the  position  of  the  drift 
object  to  obtain  the  most  accurate  tidal  current  information.  For  example,  a  4.3  knot  (215  cm/sec)  current 
will  carry  a  drifting  object  form  one  2x2  nmi  grid  square  through  a  second  and  into  a  third  in  one  hour. 
Therefore,  it  is  necessary  to  put  in  a  corrected  position  at  thirty  minute  intervals  to  insure  use  of  the  correct 
tidal  current  vector  in  forecasting  the  drift  of  an  object. 

Another  area  of  interest  in  SAR/MEP  planning  is  on  the  south  side  of  Long  Island  Sound  between  North- 
port,  73°20’W  and  Stony  Brook,  73°10’W.  In  this  area  there  occurs  a  counter  current  during  most  of  the 
flood  tide.  This  counter  current  flows  eastward  with  speeds  of  up  to  0.5  knots  (25.5  cm/s).  There  is  a  weak 
extension  of  this  counter  current  as  far  east  as  Port  Jefferson,  73’03’W  during  the  latter  stages  of  the  flood 
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tide.  A  similar  westward  flowing  counter  current  forms  during  ebb  tide  in  the  same  area.  Great  care  should 
be  taken  in  establishing  the  start  position  for  drift  forecasting  in  these  areas.  If  the  position  is  not  known 
with  certainty,  then  alternate  start  positions  should  be  considered. 
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VERIFICATION  OF  THE  LONG  ISLAND  SOUND  PROGRAM 


The  assumption  upon  which  the  LITSF  program  is  based  is  that  current  vectors  obtained  from  Roberts 
Current  Meters  at  depths  of  10  to  15  m  may  be  used  for  surface  tidal  currents.  To  test  the  validity  of 
this  assumption,  the  tidal  current  vectors  in  the  LITSF  program  for  Long  Island  Sound  are  compared  with 
surface  tidal  current  vectors  obtained  by  a  current  pole  survey  in  Long  Island  Sound  (La  Lacheur 
and  Sammons  1932).  The  current  poles  were  15  ft.  (4.6  m)  long  and  weighted  at  one  end  so  that  they 
floated  vertically  with  1  ft.  (0.3  m)  sticking  out  of  the  wate..  The  pole  is  allowed  to  drift  for  a  time 
with  a  line  attached  to  the  upper  end.  The  length  of  line  which  pays  out  is  taken  as  equal  to  the  distance  the 
pole  had  drifted.  This  is  divided  by  the  drift  time  to  determine  the  speed.  The  direction  is  determined  by 
observation  from  the  survey  ship.  This  method  gives  a  vertically  integrated  measure  of  the  upper  4.3  m 
?auer  current,  both  tidal  and  estuarine  flow. 

The  accuracy  of  the  direction  of  drift  observation  made  with  current  poles  is  estimated  to  be  ±  20°  and 
the  accuracy  in  magnitude  to  be  no  better  than  ±0.2  knots  (10  cm/sec).  The  Roberts  Current  Meters 
have  a  ±10  cm/sec  accuracy  in  magnitude  (Hicks  1967).  Agreement  between  the  two  measurements  is 
defined  as  occurring  whenever  the  error  bars  overlap.  That  is  to  say,  whenever  the  direction  is  within  30°  the 
precision  of  the  measurements  does  not  allow  us  to  identify  any  real  difference.  Similarly  whenever  the 
speeds  were  within  ±0.4  knots  (20  cm/sec)  it  was  not  possible  to  identify  any  significant  difference  in  the 
magnitude  of  the  current.  Ten  stations  were  selected  for  examination.  They  were  selected  as  areas  where 
one  might  expect  the  tidal  current  program  to  give  poor  results.  As  a  result  the  test  sample  is  not  a  random 
sample  and  should  in  fact  give  a  worst  case  evaluation  of  Long  Island  Sound.  The  surface  currents 
recorded  by  La  Lacher  and  Sammons  (1932)  are  average  currents  for  a  one  year  period.  They  are  compared 
to  the  yearly  averages  of  the  currents  in  the  Long  Island  Sound  Program. 

Station  1.  At  Throgs  Neck,  41°54’N,  73>>40,W,  the  speed  and  direction  values  were  in  agreement  100%  of 


the  time.  A  graphical  representation  of  the  magnitudes  and  directions  is  given  in  Figure  1  of  this  section. 

Station  2.  This  station  is  located  at  40°56’N,  73°27'W  is  midway  between  the  northern  and  southern 
shores  at  the  western  end  of  the  Sound.  The  speed  of  the  current  and  the  direction  of  the  L1TSF  program 
agree  with  the  test  observations  100%  of  the  time. 

Station  3.  This  station  is  located  on  the  southern  shore  at  40°58’N,  73'06’W.  This  station  is  in  the  area 
where  a  back  flow  was  predicted  in  the  LITSF  program.  The  surface  observations  also  show  evidence  of  the 
back  flow.  At  station  3,  however,  the  observed  flood  is  one  hour  longer  than  the  flood  tide  predicted  by 
the  LITSF  program,  and  the  ebb  is  one  hour  shorter.  This  produces  significant  differences  in  the  direction 
for  a  period  equal  to  1/4  of  the  tidal  cycle.  If  the  directions  are  disregarded,  then  the  magnitude  agrees  to 
the  precision  of  the  observations.  This  back  flow  appears  to  be  the  result  of  bottom  and/or  sidewall  inter- 
action  with  the  flow.  It  varies  with  depth  and  very  probably  varies  significantly  in  the  horizontal  as  well. 

Station  4.  This  station  is  located  on  the  southern  shore  at  4CP56’N,  72  '45’W.  At  this  station  there  is  also 
a  back  flow  in  the  surface,  but  it  is  different  from  that  at  station  3.  That  is,  the  observed  surface  flood 
tide  is  six  hours  long  and  the  observed  ebb  is  seven  hours  long.  The  LITSF  program  flood  is  7  hours  long 
and  the  ebb  is  six  hours  long.  As  with  the  observations  at  station  3  this  back  flow  and  temporal  variation 
is  probably  due  to  boundary  layer  effects.  The  errors  in  the  LITSF  program  could  only  be  corrected  by  a 
detailed  current  meter  project  in  the  area.  The  magnitude  of  the  currents  is  in  agreement  85%  of  the  time  at 
this  station,  and  if  the  currents  are  averaged  over  the  entire  tidal  cycle  they  agree  to  the  precision  of  the 
measurements.  The  observed  currents  are  twice  as  large  of  the  LITSF  program  currents  for  15%  of  the  time, 
one  hour  during  each  ebb  tide. 

Station  5.  This  station,  located  at  41*08^, 72°45’W  is  near  the  geographical  center  of  Long  Island  Sound. 
The  magnitude  of  the  observed  surface  currents  and  the  LITSF  program  currents  are  in  agreement  100%  of 
the  time  to  the  precision  of  the  measurements.  There  is  not,  however,  a  similar  agreement  in  the  direction. 
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The  LITSF  program  currents  are  30°  to  the  right  of  the  observed  surface  currents  during  the  flood  tide  and 
3(f  to  the  left  during  the  ebb  tide.  That  is,  the  LITSF  currents  were  northwest  and  northeast  of  the  obser¬ 
ved  surface  currents.  This  was  the  worst  case  for  directional  agreement  and  the  only  one  with  complete 
disagreement. 

Station  6.  The  observed  surface  currents  and  the  LITSF  currents  at  this  station,  41°14’N,  72°42’W, 
agreed  in  both  magnitude  and  direction  100%  of  the  time. 

Station  7.  This  station  is  off  the  mouth  of  the  Connecticut  River  at  41°16’N,  72°21’W.  The  directions 
agree  within  the  precision  of  the  measurements.  The  direct  observations  of  the  surface  currents,  however, 
show  that  on  the  surface  the  current  begins  to  ebb  about  one  hour  before  it  does  at  the  mid-level  depth  of 
the  current  meter  used  to  generate  the  tidal  currents  for  the  LITSF  program.  Such  vertical  variations  are 
not  uncommon  in  river  channels  and  at  the  mouths  of  rivers.  The  directions  agree  to  the  precision  of  the 
measurement,  except  for  the  last  hour  of  flood. 

Station  8.  This  station  is  located  midway  between  the  northern  and  southern  shores  at  the  eastern  end  of 
Long  Island  Sound  at  41°04’N,  72°15’W.  The  observed  surface  currents  and  the  LITSF  program  currents 
agree  in  magnitude  and  direction  100%  of  the  time. 

Station  9.  This  station  is  on  the  Connecticut  shore  just  west  of  the  mouth  of  the  Themes  River 
at  41°  16’N,  72?>06’W.  The  LITSF  program  currents  agree  in  direction  85%  of  the  tidal  cycle.  The  magnitude 
agrees  with  surface  observations  throughout  the  flood  tide,  but  on  the  ebb  tide  the  LITSF  program  eurrents 
are  180%  of  the  observed  values.  This  overestimation  lasts  for  three  hours  or  approximately  23%  of  the  total 
tidal  cycle.  As  in  the  case  of  station  7  these  changes  are  primarily  due  to  the  effects  of  the  river  outflow. 

Station  10.  Station  10  is  at  the  Race,  41°14’N,  72?>03’W.  The  LITSF  program  currents  agreed  in  direction 
92%  of  the  tidal  cycle  when  compared  with  surface  observations.  The  LITSF  program  currents  were  faster. 
110%,  than  the  observed  values  during  the  flood.  During  the  ebb  currents,  the  LITSF  program  currents 
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were  90%  of  the  observed  surface  currents.  The  differences  amounted  to  a  half  a  knot  in  absolute  values 


on  the  average. 


DISCUSSION 


Although  this  test  was  designed  to  give  a  worst  case  view  of  the  LITSF  program,  it  showed  that  complete 
agreement  in  speed  and  direction  was  observed  83.2%  of  the  time.  At  40%  of  the  stations  the  agreement  was 
100%  of  the  time.  Only  one  station  was  off  for  the  entire  tidal  cycle.  That  station  was  3(f  off  in  direction 
for  the  entire  tidal  cycle.  Disagreement  appears  to  occur  most  often  near  the  coast  and/or  near  river  mouths. 
In  both  cases,  the  influence  of  the  bottom  topography  and  variable  river  flow  are  large.  The  bottom  topo¬ 
graphic  features  such  as  sand  bars  are  also  subject  to  a  relatively  slow  variation  with  time.  This  makes 
correcting  the  model  for  such  influences  a  very  difficult  process  and  one  of  limited  value.  Because  both  the 
LITSF  values  and  the  current  pole  survey  are  direct  current  observations,  the  differences  are  to  some  extent 
a  measure  of  the  natural  variability  in  the  currents. 

CONCLUSIONS 

The  use  of  mid  level  current  observations  to  infer  surface  currents  is  apparently  a  valid  assumption.  It 
is  recommended  that  the  currents  from  the  National  Ocean  Survey  Tidal  Current  Charts  for  Long  Island 
Sound  be  used. 
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SPECIAL  FEATURES  OF  SAN  FRANCISCO  BAY 

Although  San  Francisco  Bay  is  an  extensive  embayment,  it  possesses  only  one  entrance  located  at  the 
Golden  Gate  region  (Figure  2).  In  addition  to  this  restricted  entrance,  the  bay  can  be  divided  into  two 
distinct  tidal  regimes;  a  standing  and  a  progressive  tidal  wave.  The  Golden  Gate  and  Southern  Bay  region, 
contain  a  tidal  scheme  similar  to  a  standing  wave,  as  evidenced  by  simultaneous  changes  in  the  tide  through¬ 
out  the  area  (Disney  et.  al.  1925).  The  northern  bay  region  and  San  Pablo  Bay  exhibit  a  progressive  wave 
phenomenon.  This  occurance  shows  up  as  a  progresively  later  flood  or  ebb  current  as  one  moves  north  to 
San  Pablo  Bay.  The  time  of  occurance  of  these  currents  also  lags  the  southern  regions’  current  changes. 

In  the  San  Francisco  Bay  system,  friction  is  the  most  significant  factor  modifying  the  tidal  currents. 
This  shows  up  vividly  in  the  Golden  Gate  area  (bounded  by  (1)  37°  48’N  -  122°  22.5"W  (2)  37°50’N 
122 °32.5'Vil  (3)  37°48’N  122°-  32.3^  and  (4)  37°50’N  -  122°  22.5’W)  as  a  very  pronounced  horizontal 
current  variability.  This  variability  presents  itself  as  a  set  of  eddies  formed  by  lateral  friction  on  the  northern 
and  southern  shores  near  the  western  approach  to  Golden  Gate.  These  eddies  produce  near  shore  counter 
currents  on  both  the  ebb  and  flood  tides.  As  a  result,  a  drifting  object  near  shore  and  just  seaward  of  the 
Golden  Gate  Bridge  may  travel  westward  on  a  flood  tide  or  into  the  bay  on  an  ebb  tide.  This  feature  makes 
it  extremely  important  to  know  the  initial  position  and  start  time  for  any  SAR  case  when  predicting  a  drift 
path. 

The  duration  of  flood  and  ebb  tides  as  well  as  the  magnitude  of  the  currents  is  also  found  to  be  variable 
within  the  bay.  The  flood  to  ebb  duration  was  found  by  Disney  (1925)  to  be  6.71  hours  to  5.71  hours 
(1.18:1).  This  ratio  persists  throughout  most  of  San  Francisco  Bay  with  the  exception  of  San  Pablo  Bay 
where  the  influence  of  river  run  off  causes  a  revised  ratio.  Also,  in  the  northern  portion  of  the  Bay  and  San 
Pablo  Bay,  the  initial  flood  current  is  sub-surface  wand  moves  vertically  to  the  surface,  while  the  ebb  exhibits 
the  reverse  (Disney  1925).  The  magnitude  of  the  currents  through  the  Bay  at  times  exceeds  2  knots 
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(102.1  cm/sec).  In  an  area  of  great  horizontal  variability,  they  may  be  as  high  as  4  knots  (204  cm/sec. 

To  model  the  San  Francisco  Bay  current  fields  with  the  above  noted  variability,  a  flexiable  grid 
pattern  and  set  time  interval  were  selected.  In  areas  of  high  current  variability,  a  small  grid  square  (lxl 
nautical  miles)  was  utilized,  while  a  larger  grid  (2x2  nautical  miles)  was  set  up  in  the  rest  of  the 
Bay  (Figure  3).  Even  with  the  noted  inequality  of  the  ebb  to  flood  time,  the  tidal  sequence  was  divided  into 
equal  time  intervals.  These  modeled  tidal  hours  are  equivalent  to  solar  hours  set  at  1.12  hours  for  the  flood 
and  0.95  hours  for  ebb.  These  parameters  of  grid  size  and  time  interval,  do  however,  require  that  in  drift 
studies  the  calculated  tidal  current  vector  be  updated  on  an  hourly  basis  or  less  if  in  an  area  of  high  currents. 
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VERIFICATION  OF  THE  SAN  FRANCISCO  BAY  PROGRAM 
To  establish  the  validity  of  the  total  current  vectors  calculated  by  the  San  Francisco  Bay  computer 
program,  comparisons  between  calculated  and  measured  values  were  made  at  fifty  one  stations  in  tne  Bay 
(Figure  3a  and  3b).  The  measured  current  values  were  obtained  from  the  hydraulic  tests  run  on  the  l!.S. 
Army  Corps  of  Engineers,  San  Francisco  Bay  Delta  Model,  Sausulito,  California.  These  tests  were  run  in  May 
1978  during  a  simulation  of  the  1977  Delta  dynamics  involving  a  net  river  flow  (Sacramento  and 
San  Joaquin  Rivers)  of  4700  cubic  feet  per  second  (133.09  m  /sec),  an  ocean  salinity  of  33.00%o  and  a  19 
year  mean  tidal  curve  as  the  repetitive  forcing  function  near  the  mouth  of  San  Francisco  Bay.  The  com¬ 
puterized  tidal  current  vectors  utilized  in  the  verification  were  calculated  for  three,  twenty  four  hour  time 
periods  at  the  fifty  one  sampling  stations.  The  time  periods  were  selected  so  that  the  computer  generated 
currents  would  be  in  two  catagories:  (1)  a  set  in  which  the  tidal  height  curve  near  the  mouth  of  the  Bay 
approached  the  19  year  mean  and  (2)  a  set  in  which  the  tidal  current  curve  near  the  mouth  of  the  Bay 
approached  the  Delta  Model  Tidal  Current  values  measured  at  the  Golden  Gate  Bridge.  The  first  set  was 
obtained  from  two  time  periods;  10-11  March  and  14-15  May  1977  (Figures  4  and  5)  and  the  second  was 
matched  with  a  period  during  16-17  October  1977  (Figure  6). 

COMPARATIVE  METHOD 

The  degree  of  agreement  between  the  observed  and  calculated  tidal  currents  at  the  sampling  stations  was 
analyzed  using  portions  of  the  Student-t  formula  (equation  1,  Peterson,  1973)  and  an  hourly  comparison 
of  the  speed  and  direction  values.  Normally,  the  calculation  of  the  Student-t  values  for  a  paired  data  set 
would  involve  the  construction  of  a  hypothesis  test  to  determine  the  probability  that  the  sample  means 
differed  by  some  set  value  (in  this  study  the  value  was  set  at  zero  in  equation  1).  However,  this  course  of 
evaluation  would  have  been  inconclusive,  because  even  through  the  difference  in  the  means  approaches 
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Figure  3a.  Northern  Section  of  San  Francisco  Bay  showing  partial  grid  layout  of  current  file  and  sampling 
stations  1-37. 
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Figure  3b.  Southern  Section  of  San  Francisco  Bay  showing  partial  grid  layout  of  current  file  and  sampling 
stations  38-51. 


SAN  FRANCISCO  BAY- DELTA 
MODELED  CURVE 


O 


Figure  5.  Modeled  and  predicted  tidal  height  curve  at  San  Francisco’s  Golden  Gate.  Modeled  curve 
represents  19  years  mean  and  predicted  curve  for  period  1500  14  May  to  1600  15  May  1977. 


PREDICTED  TIDAL 
CURRENT  CURVE 


TIME  IN  HOURS 

Modeled  and  predicted  tidal  current  curve  at  San  Francisco’s  Golden  Gate.  Modeled  curve 
represents  measured  values  from  Delta  Bay  Hydroulic  Model  1977  Delta  experiment  and 
predicted  curve  for  period  0624  16  October  to  0724  17  October  1977. 
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The  standard  deviation  (B)  was  then  used  to  calculate  the  standard  error: 

c=jTT 

where: 

f -  1 

I  1  +— the  adjustment  value  required  when  using  the  standard  deviation  in  the  calculations. 

V  n|  nX 

(B)  =  the  standard  deviation  from  equation  1 

(C)  =  the  standard  error. 

The  non  statistical  methodology  utilized  in  the  hourly  comparisons,  consisted  of  determining  the  number 
of  times  per  complete  tidal  cycle  that  the  computer  model's  speed  or  direction  values  were  higher/lower 
or  left/right  of  the  hydraulic  model  values.  These  values  were  then  computed  as  a  percentage  and  a  mean 
value  for  all  catagories. 

Once  all  station  values  both  statistical  and  general,  were  calculated,  these  were  grouped  into  four  geo¬ 
graphical  areas  (Table  la;  and  Figure  2)  and  a  quantitative  set  of  values  was  determined  for  the  region  (Table 
lb).  The  regional  set  of  values  consisted  of:  (1)  a  mean  difference  range;  (2)  the  standard  error  range: 
(3)  percentage  range  of  computer  values,  high  or  low,  left  or  right  of  the  hydraulic  model;  (4)  mean  per¬ 
centage  of  occurrences  high/low,  or  left/right  of  the  hydraulic  model  and  (5)  mean  value  for  occurrance 
high/low  or  left/right  of  hydraulic  model. 
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TABLE  la 

Geographical  grouping  of  sampling  stations  in  San  Francisco  Bay  for  analysis 


STATIONS 


BOUNDARIES  GENERAL  DESCRIPTION 


01-09 


37.78°  to  37.8?N 
122.50°  to  122.33^ 


Entrance  to  San  Francisco 
Bay 


10-19 


37.80°  to  37.8?  N 
122.33°  to  122.42°W 


Area  surrounding  Yerba  Beuna 
Islands 


TABLE  lb 


n 


zero,  the  shapes  of  the  two  curves  could  be  completely  different.  Therefore,  the  Student’s-t  formula 
(equation  1)  was  used  only  to  determine  a  mean  difference  and  a  standard  error  for  each  set  of  paired  station 
data.  The  above  noted  values  for  each  set  of  paired  station  data  were  determined  using  the  Student’s-t 
statistics  evaluation  program  (comparison  of  population  means)  as  contained  in  the  library  programs  of 
Texas  Instruments  TI  58/59  programmable  calculator  (Applied  Statistics,  1977).  The  program  involved  the 
use  of  the  bivariate  data  entry  program  and  a  two  sample  test  to  calculate  the  individualized  student-t  values 
of  each  station,  the  difference  of  the  means  and  then  the  standard  deviation,  which  was  adjusted  to  be  the 
standard  error.  The  equations  utilized  by  this  TI  58/59  program,  as  stated  by  Perterson  (1973)  are: 


where: 

Z 

X 

X 

M 

N, 

Na. 
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Student-t  value  with  nj  +  n^-2  degrees  of  freedom 

24  hour  means  for  the  computer  generated  currents/direction 

24  hour  mean  for  the  hydraulic  model  currents/direction 

hypothetical  mean  difference  value;  set  at  zero 

number  of  computer  model  observations 

number  of  hydraulic  model  observations 

and  n(  =  n^ 

standard  deviation 

difference  between  the  means 


J 


(1) 
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STATIONS 

BOUNDARIES 

GENERAL  DESCRIPTION 

20-37 

37.83°  to 

38.06°N 

Northeastern  San  Francisco 

122.48°  to 

122.29^ 

Bay  and  San  Pablo  Bay 

38-52 

37.50°  to 

37.70°N 

Southern  San  Francisco  Bay 

122.35°  to 

122.13^ 

RESULTS 

The  statistical  and  general  analysis  values  (Table  2  )  show  that  the  two  models  differ  significantly  in  the 
magnitude  and  direction  of  the  tidal  current  vectors  each  produce.  The  range  in  the  speed  and  direction 
mean  differences  were  found  to  be  approximately  equal  for  the  months  of  March  and  May  for  speed  and 
all  three  months  for  direction,  but  they  differed  in  magnitude  region  by  region.  The  observed  speed  ranges 
were  (largest  to  smallest  range): 


TABLE  2 

Statistical  and  general  analysis  values  for  the  four  regions  of  Figure  2 
REGION  MONTHS  VALUE  RANGE  IN  KNOTS 


I 

All 

-.82  to  .87 

IH 

March  and  May 

-.28  to  .37 

All 

-.33  to  .41 

n 

March  and  May 

-.15  to  .51 

All 

-.65  to  .64 

IV 

March  and  May 

-.15  to  .31 

All 

-.76  to  .32 

The  relative  size  of  the  overall  range  in  the  mean  direction  differences  did  not  coincide  with  the  speed 
distribution.  Rather,  the  distributions  from  the  largest  to  the  smallest  were: 
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TABLE  3 

Variations  in  direction  in  the  four  regions  of  Figure  2 


REGION  MONTHS  VALUE  RANGE  IN  DEGREES 


III 

All 

-32.78 

to 

110.03 

IV 

March  and  May 

-61.29 

to 

13.84 

All 

-61.29 

to 

47.30 

II 

All 

-74.11 

to 

11.94 

I 

All 

-56.50 

to 

-2.70 

To  further  expand  how  significant  the  mean  difference  values  are,  an  examination  of  the  overall  range  in 
the  standard  error  is  necessary.  This  gives  some  indication  as  to  the  agreement  between  a  set  of  curve  shapes 
and  to  the  overall  agreement  between  the  two  models.  Examining  the  regions  for  speed,  the  greatest  to  the 
least  range  in  standard  error  are: 

TABLE  4 

Variations  of  magnitude  in  the  four  regions  of  Figure  2 


REGION 

RANGE  IN 

KNOTS 

II 

0.06  to 

1.67 

IV 

0.14  to 

1.02 

I 

0.10  to 

0.95 

HI 

0.09  to 

0.83 

The  sequence  of  the  standard  error  range  for  direction  varied  from  the  speed  distributions  as  follows: 


TABLE  5 

The  standard  error  range  for  direction  in  the  four  regions  of  Figure  2 


REGION 

I 

m 

IV 

n 


RANGE  IN  DEGREES 

17.60  to  45.20 

19.55  to  46.69 

6.73  to  28.79 
25 

23.76  to  29.69 


Even  though  these  values  appear  to  be  small  in  magnitude,  overall  0.09  to  1.02  knots  (4.6  to  52.0 
cm/sec)  for  speed  and  6.73°  to  46.69°  for  direction,  they  indicate  that  a  significant  disagreement  in  curve 
shapes  or  values  exists.  To  more  clearly  determine  how  large  the  disagreement  is  between  the  computerized 
and  the  hydraulic  model,  the  curve  values  were  compared  on  an  hour  by  hour  basis.  The  results  of  Tables  4 
and  5;  are  further  summerized  below  as  an  overall  percentage  of  occurrence  and  a  mean  value  for  a  specific 
region  to  give  a  clearer  picture  of  agreement/disagreement.  These  values  were  found  to  be: 

TABLE  6 

Percentage  of  high  and  low  magnitude  values  relative  to  the  mean 
SPEED  (KNOTS) 


REGION 

HIGH  %  OF  OCCURRANCE 

VALUE 

LOW  %  OF  OCCURRANCE 

VALUE 

I 

48.33 

0.91 

51.67 

-1.01 

II 

40.30 

0.38 

59.70 

-0.56 

III 

55.00 

0.60 

45.00 

-0.71 

IV 

48.33 

0.50 

51.67 

-0.54 

Percentage  of  left  and  right  occurence  relative  to  the  mean 
DIRECTION 

REGION 

RIGHT  %  OF  OCCURRANCE 

VALUE 

LEFT  %  OF  OCCURRANCE 

VALUE 

I 

77.67 

62.23 

22.33 

-63.11 

II 

76.33 

66.11 

23.67 

-71.05 

m 

75.33 

49.15 

24.67 

-75.81 

IV 

89.33 

46.89 

10.67 

-71.21 

It  can  be  seen  that  speed  is  under  predicted  in  Regions  I,  II,  and  IV  over  50%  of  the  time  by  0.5  to  1.0 
knots  (25.5  to  51.479  cm/sec)  and  over  predicted  in  Region  III  over  50%  of  the  time  by  0.6  knots  (30.6 
cm/sec).  Where  as,  direction  is  predicted  high  or  to  the  right  of  the  observed  values  over  75%  of  the  time  by 
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values  ranging  from  49.15°  to  66.11°.  The  speed  curves  show  that  the  absolute  value  of  the  current 
magnitudes  varied  from  .06  knots  to  3.38  knots  (3.1  to  172.4  cm/sec).  Thus  the  predicted  difference 
was  high  18  to  100%  depending  on  the  magnitude  of  current  vector.  Since  direction  is  not  an  increasing 
magnitude,  the  sector  of  agreement/disagreement  ranged  from  118.10°  to  137.1^  in  width,  which  far  out 
ways  any  usefulness. 
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DISCUSSION 


The  development  of  the  computerized  San  Francisco  Bay  tidal  current  model  was  an  attempt  to  create  a 
quick  access  reference  of  the  tidal  current  at  any  point  or  time  in  the  Bay  for  application  to  Search  and 
Rescue  drift  programs.  The  statistical  and  general  comparisons  of  the  tidal  data  produced  by  the  model 
to  the  hydraulic  model  of  the  Army  Corps  of  Engineers  was  an  attempt  to  determine  if  the  two  data  bases 
were  comparable.  The  Bay  hydraulic  model  was  utilized  as  a  reference  source  to  simplify  data  collection 
(no  ship  time  or  current  meters  required)  and  analysis  by  eliminating  any  wind  drift  vectors  produced  in 
real  world  sampling. 

The  San  Francisco  Bay  Hydraulic  Model  is  driven  by  a  repetitive  19  year  mean  tidal  forcing  function  at 
the  ocean  entrance  and  comes  to  equilibrium  throughout  the  Bay  after  spin-up.  This  allows  for  repetitive 
sampling  of  data  curves,  but  does  not  match  real  world  conditions  since  the  real  ocean  entrance  tidal  curve 
is  continually  being  modified  by  tidal  constituents.  To  over  come  part  of  this  problem,  the  ocean  entrance 
curves  were  matched  for  tidal  height  (March  and  May)  and  tidal  current  speeds  (October)  and  the  assumption 
made  that  all  other  geographic  locations  may  approach  agreement.  However,  as  seen  by  Table  lb,  this 
was  not  the  case.  This  can  be  seen  in  the  large  speed  and  direction  disagreements.  Another  reason  for  this 
nonagreement  is  felt  to  be  due  to  the  fact  that  the  computerized  tidal  current  charts  are  based  on  a  1954 
tidal  current  survey  by  the  National  Ocean  Survey.  Therefore,  comparisons  of  a  1954  survey  to  1977  tidal 
conditions,  as  modeled  by  the  hydraulic  model,  should  show  disagreement.  NOS  is  in  the  process  of  resur¬ 
veying  the  Bay  and  plans  for  completion  of  field  work  in  1980. 
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CONCLUSIONS 


The  comparison  between  the  data  sets  of  the  computerized  tidal  current  model  and  the  hydraulic  San 
Francisco  Bay  Model  revealed  that  a  significant  disagreement  in  speed  and  direction  exists.  It  is  recom¬ 
mended  that  the  National  Ocean  Survey  tidal  current  chart  data  be  used  with  caution.  The  disagreement 
in  the  tidal  current  speed  was  found  to  be  best  represented  by  the  hour  by  hour  analysis,  which  showed 
the  computer  model  underpredicting  speeds  and  overpredicting  direction.  The  speed  values  were  under- 
predicted  in  magnitude  in  Regions  I,  II  and  IV  over  50%  of  the  time  by  from  0.5  to  1.0  knots  (25.5  to 
51.0  cin/sec)  and  overpredicted  in  Region  III  over  50%  of  the  time  by  0.6  knots  (30.6  cm/sec).  The  under- 
prediction  of  speed  was  found  to  be  approximately  18%  to  1000%  of  the  magnitude  of  a  current  generated 
by  the  San  Francisco  Bay  hydraulic  model,  another  indication  of  significant  disagreement.  The  direction 
values  were  predicted  high  or  to  the  right  of  an  observed  value  over  75%  of  the  time  for  all  regions  by  a  value 
between  49.15°  to  66.11°  .  The  Coast  Guard  Oceanographic  Unit  is  planning  a  field  test  of  the  accuracy 
and  utility  of  the  Army  Corp  of  Engineers  Hydraulic  Model  and  the  National  Ocean  Survey  Tidal  Current 
Charts  as  sources  of  tidal  current  information.  Should  this  field  test  show  the  Hydraulic  Model  to  be  a 
superior  source  of  tidal  current  data  a  new  set  of  tidal  current  charts  will  be  generated  based  on  it. 
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APPENDIX  I 


INSTRUCTIONS  FOR  RUNNING  THE  PROGRAM 

The  program  LITSF  and  the  current  velocity  by  position  grid  are  generally  stored  on  disk.  The  following 
program  is  required  to  obtain  the  tidal  current  information  from  the  tidal  current  program  for  Long  Island 
Sound  (See  listing  at  the  end  of  this  section). 

The  first  six  cards  are  called  the  job  control  cards,  utilized  by  the  Coast  Guard’s  CDC  3300  computer. 
They  will  remain  unchanged  for  every  job  run. 

Card  LITS  I  is  the  station  name  and  meridian  card.  It  is  always  the  same  for  Long  Island  Sound. 

Card  LITS  2  contains  the  permanent  current  indexing  parameters  and  directions  for  Long  Island  Sound. 
These  are  entered  in  the  following  format:  Columns.  1-6  Permanent  current  (PERMC)  3  decimals.  Col¬ 
umns  7-8  IND1,  Columns  9-10  IND2,  Columns  11-12  IND3,  Columns  13-14  IND4,  Columns  15-16  IND5, 
Columns  17-18  IND6,  Columns  19-22  Flood  Direction  (NFDIR),  Columns  23-26  Ebb  Direction  (NEDIR). 
The  function  of  each  of  the  indexing  parameters  is  given  in  Table  1  at  the  end  of  this  section.  Normally 
these  indexing  parameters  are  never  changed. 

Cards  LITS  3-8  are  the  station  constants  which  determine  the  amplitude  and  phase  lag  for  the  reference 
station.  These  constants  are  supplied  by  the  National  Ocean  Survey  Production  Division  and  need  not  be 
changed.  The  constitutents  are  entered  in  the  following  format:  Columns  1-4  station  number,  Columns  5-8 
station  card  sequence  number,  Colums.  9-13  constituent  amplitude  (3  decimal  places),  Columns.  14-17 
constituent  epoch  (l  decimal  place),  Columns  18-71  repitition  of  Columns  9-17. 

Cards  LITS  9-23  are  the  year  constant  cards.  The  format  is  given  by  Morgan,  et  al.  1975  is  as  follow: 
“Columns  1-4  year-example  =  1979,  Columns  5-6  year  type  identifer  =  0  for  a  common  year,  =  1  for  a  leap 
year,  Columns  7-8  year  card  sequence  number,  Columns  9-12  Node  factor  (f)  for  the  first  constituent, 
Columns  13-16  Greenwich  V  +  u  for  the  first  constituent  (Oh,  Jan.  1),  Columns  17-72  repititions  of  columns 
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9-16,  Each  year  the  node  factor  (f)  and  equilibrium  argument  (V  +  u)  must  be  added  for  the  following 
twenty  six  constants:  M2  ,S  2,N2  ,K|  ,M*  ,0)  ,V1« ,(MK)3,S4,(MN)4,  ^.S  6,P|(2!N)2,(00)i  ,^?Si  ,Mi,Ji,Mm, 
Ssa.Sa,MSf,Mf.  p  ,Q ,  ,T  2  ,Rj  ,(2Q)  1  ,P  i,(2SM)a,M3,L  2,(2Mk)3,Mg  and(MS)4.”  These  constituents  are 
found  in  (Schureman  1958). 

f.ard  1.1TS  24  is  the  Date  Control  Card  which  is  the  first  card  which  the  user  must  supply.  This  card 
is  to  have  the  following  format:  Columns  1-2  month  number,  Columns  3-4  beginning  day,  Colums  5-6 
ending  day.  Columns  7-72  repetition  of  the  above  as  needed.  The  maximum  period  for  a  single  calculation 
is  32  days.  If  the  desired  time  periods  are  not  consecutive,  use  individual  date  control  cards  in  separate 
computer  runs. 

Card  25  is  the  Termination  Control  Card.  This  card  determines  if  the  tides  will  be  calculated  for  one  or 
more  than  one  reference  station  or  time.  It  has  the  following  format:  Columns  1-4  (MS)  =  0.  The  new 
problem  is  for  the  same  station;  =  1.  The  new  problem  is  for  a  different  station.  Columns  5-8  (MY)  =  0.  The 
new  problem  is  for  the  same  year;  =  1.  The  new  problem  is  for  a  different  year,  Columns  9-12  (MD)  =  0. 
The  new  problem  is  for  the  same  day;=  1.  The  new  problem  is  for  a  different  day.  The  job  is  terminated  if 
all  three  cards  are  zero.  This  is  the  format  which  is  always  used  in  calculating  tidal  currents. 

Card  LITS  26  this  card  may  be  used  to  call  for  a  listing  of  the  tidal  currents  at  the  reference  station  for 
the  period  given  in  the  Date  Control  Card.  If  this  is  desired  a  1  is  entered  in  column  1.  If  not  a  0  is  entered 
in  column  1.  When  using  the  model  to  calculate  tidal  currents  a  0  is  entered  in  the  column. 

Cards  26-29  are  repeated  for  each  additional  time  or  position  for  which  a  tidal  vector  is  needed. 

Cards  UTS  30  is  the  last  Month  List  card  it  is  always  blank. 

Card  UTS  31  is  the  final  Day  Count  and  Station  I.D.  card. 
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Card  LITS  32  is  the  last  position  card.  A  latitude  of  99  is  entered  in  columns  one  and  two.  This  signals 
the  end  of  the  input  data  cards. 

The  last  two  cards  in  the  deck  are  the  job  termination  control  cards  used  in  the  Coast  Guard’s  CDC 
3300,  they  never  change. 
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APPENDIX  II 


INSTRUCTIONS  FOR  RUNNING  THE  MODEL  FOR  SAN  FRANCISCO  BAY 

The  following  small  program  is  required  to  run  the  tidal  current  model  for  San  Francisco  Bay.  The  listing 
is  identical  to  the  Long  Island  Sound  program  for  all  job  control  cards,  but  exhibits  necessary  changes  for 
program  data  and  instructions. 

The  first  six  cards  are  called  the  job  control  cards.  They  will  remain  constant  for  every  job  run. 

Card  1.  Selects  the  reference  tidal  station  and  the  time  meridian  longitude.  It  is  always  the  same  for 
San  Francisco  Bay. 

Card  2.  Contains  the  permanent  current  indexing  parameters  and  directions  for  San  Francisco  Bay.  The 
format  used  is  the  same  as  that  used  in  the  program  for  Long  Island  Sound  (See  Instructions  For  Running 
The  Program  For  Long  Island  Sound  Section). 

Card  3-9.  Are  the  station  constants  consisting  of  the  amplitude  and  phase  lag  for  each  tidal  forcing 
function  at  the  reference  station.  The  format  is  the  same  as  that  in  the  program  for  Long  Island  Sound. 

Cards  9-23.  Are  the  year  constant  cards.  The  format  for  these  cards  and  method  of  determining  these 
constituents  is  outlined  in  the  section  giving  instructions  for  running  the  model  for  Long  Island  Sound. 

Card  24.  Is  the  date  control  card.  Its  purpose  is  to  set  the  month(s),  start  and  stop  day  that  tidal  currents 
are  to  be  computed  for.  It  has  the  same  format  as  Card  24  for  Long  Island  Sound  (See  Appendix  I). 

Card  25.  Is  the  termination  control  card  which  is  used  to  separate  different  reference  stations  and  time 
controlled  calculations.  The  format  for  this  card  is  given  in  the  section  giving  instructions  for  the  running 
of  the  program  for  Long  Island  Sound. 

Card  26.  Is  the  month  list  card.  This  card  specifies  that  a  listing  of  the  times  of  slack  water  maximum 
ebb  and  flood  for  the  interval  given  on  the  date  control  card  is  to  be  printed.  A  one  in  column  (I)  causes 
only  a  listing  to  be  printed.  A  zero  causes  the  program  to  calculate  the  tidal  current  at  the  position  of 
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Cards  25  and  26.  Are  used  to  consecutively  run  two  or  more  differing  tidal  current  information  request 
for  information  at  reference  stations.  They  are  normally  blank  when  being  used  in  operational  work. 

Card  27.  Is  the  station  ID  and  Day  count  card.  The  format  and  use  of  this  card  is  described  in  Appen¬ 
dix  I. 

Card  28.  Is  the  position  card  specifying  the  geographic  position  where  tidal  current  information  is 
required.  The  latitude  is  given  to  the  nearest  minute  and  the  longitude  to  the  nearest  0.25  minute.  The 
format  is  given  in  the  section  on  Long  Island  Sound. 

Card  29.  Is  the  date  time  card.  This  card  inputs  the  date  time  group  of  the  occurrence.  Its  format  is 
given  in  the  section  on  Long  Island  Sound. 

The  set  of  cards,  26-29.  Are  repeated  for  each  additional  time  or  position  for  which  a  tidal  current 
vector  is  needed. 

Card  30.  Is  the  last  month  list  card  and  is  placed  at  the  end  of  all  series  of  time  and  position  cards.  It  is 
always  blank. 

Card  31.  Is  the  last  day  count  and  station  ID  card.  It  is  formated  the  same  as  card  27. 

Card  32.  Is  the  final  position  card.  A  latitude  of  99  degrees  is  always  entered  in  columns  1  and  2.  This 
signals  the  end  of  the  data  end  put  cards. 

The  final  two  cards  are  the  job  termination  control  cards,  and  are  formated  as  shown.  These  two 
cards  like  the  first  six  cards  are  always  the  same. 
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APPENDIX  III 


PRODUCTION  OF  THE  TIDAL  CURRENT  FILES 

The  tidal  current  files  are  prepared  using  tidal  current  charts  form  the  National  Ocean  Survey  (NOS). 
The  most  recent  charts  should  be  used  in  making  these  files.  The  charts  are  available  from  the  National 
Ocean  Survey  Office,  Rockville,  Maryland.  The  first  step  is  the  construction  of  a  grid  overlay  for  the  tidal 
current  chart.  These  are  either  done  using  a  template  or  drawn  in  with  pencil  so  that  the  lines  lie  on  whole 
minutes  of  latitude  and  longitude. 

A  grid  square  measuring  approximately  2x2  nautical  miles  has  been  found  to  work  best.  This  grid  size 
corresponds  to  a  2’  x  3’  along  much  of  the  U.S.  coast  line.  A  tidal  current  speed  and  direction  is  assigned 
to  each  grid  square  for  the  entire  chart  of  currents.  These  values  are  either  based  on  direction  measure¬ 
ments  or  from  interpolation.  When  an  interpolated  value  is  required,  the  speed  and  direction  of  the  currents 
are  assumed  to  change  linearly  with  changes  of  equal  value  between  each  grid  square.  That  is,  if  there  are 
5  blank  grid  squares  between  two  observations,  then  the  change  from  one  block  to  the  next  will  be  1/5 
of  the  difference  between  the  two  observations.  The  same  is  true  for  the  directions.  An  example  is  given 
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The  current  is  always  assumed  to  parallel  the  coast  line  unless  observed  to  do  otherwise.  Upon  completion 
of  each  tidal  chart  the  current  values  are  transferred  from  the  array  set  up  to  a  computer  listing  with  one 
computer  card  for  each  grid  square.  The  cards  have  the  following  format: 

Line  1-4.  The  latitude  in  degrees  and  minutes  for  the  south  east  corner  of  the  grid  square. 

Line  6-9.  The  longitude  in  degrees  and  minutes  for  the  south  east  corner  of  the  grid  square. 

Line  11-14,  The  speed  in  10th  of  knots  and  direction  in  10s  of  degrees  for  slack  water  at  the  reference 
station. 


Line  16-19.  The  same  as  lines  11-14  for  one  hour  after  slack  water  at  the  reference  station. 

This  sequence  may  be  repeated  for  each  interval  of  the  tide  up  to  13  times.  All  multiples  of  5  are  to  be 
left  as  blanks. 

A  program  is  produced  by  copying  LISND  and  changing  only  the  five  character  name  to  the  first 
five  characters  in  the  name  of  the  body  of  water  being  modeled. 


fjnq.r,C0"P6??>  199!  7«SE»  1  .=-0 
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i*DEE ( A. .WKL YTMPF .LISNO, 01 .CGOu.OUOl . 80, 1000*. «  *84 1.16) 

S*DEF (0,  *  1 ,wKL YTMPF ♦LISND.01 ,C60U»0) 

AP  =  N 

SFTNU(X) 

PROGRAM  LISNO  LSNU 

DIMENSION  CUR(IOOO) .MOIM(IOOO) 

INTEGER  CUR 
INTEGER  mo  IP 
N  =  0 

SO  READING. 3)  LAf)«LAM«LOD  «  LOm  ,  ( CUR  (  I)  *  MD I K  (  I )  *1  =  1  .13)  »  I  AR 
WRITE (1.3)  LAD.LAM.LOD.LOM, (CUP(I) . MD I R ( I ) , I  =  1 . 1 3 ) . I  AD 
3  FORMAT  (  I?,  12, 13,12, 1 X, 13 < 12, I?. 1 <)  ,  I?) 


N  =  N+  1 

IF  (LOP  .FO.  99R)  GO  TO  POO 
GO  TO  M) 

200  WRITE  (61,260) 

260  FORMAT  (1h  ,*DATA  FOP  LONG  ISLAND 
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350  FORMAT  (16  .16,*  CAkDS  READ*) 
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F  N  0 

FI*' I'. 

S >  t  L  GO 
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3749 

12225  309 

1309 

1609 

1509 

1109 

409 

1027 

2127 

2427 

22?  7 

1727 

727 

1 

3749 

122262500226 

1508 

2608 

2908 

2608 

1708 

526 

2126 

3226 

3726 

3126 

1826 

1 

3750 

12225  7?ft 

506 

908 

808 

608 

30  8 

533 

1726 

2326 

2626 

2326 

1726 

1 

3750 

122262500323 

505 

905 

805 

605 

305 

533 

1623 

2323 

2323 

2023 

1423 

1 

3781 

12225  227 

309 

409 

409 

309 

109 

227 

427 

627 

727 

627 

327 

1 

3751 

122262500103 

1203 

1903 

2203 

2003 

1303 

121 

1021 

1721 

2021 

1721 

921 

1 

375? 

12225  518 

1004 

1204 

1204 

904 

404 

418 

1418 

2118 

2518 

2318 

1818 

1 

3754 

12225  618 

600 

1400 

1800 

1700 

800 

800 

218 

1218 

1818 

1718 

1318 

1 

375ft 

12225  718 

200 

1200 

1800 

2000 

1600 

900 

218 

1518 

1818 

2018 

1718 

1 

3758 

12225  2019 

103 

1203 

2103 

2303 

2003 

1203 

219 

1919 

3119 

3419 

3010 

1 

3800 

12225  620 

502 

110? 

1402 

160? 

1502 

80? 

320 

1520 

2020 

2020 

1820 

1 

3802 

12225  620 

502 

1102 

1402 

1602 

150? 

80? 

320 

1520 

2020 

2020 

1820 

1 

3804 

12225  315 

533 

833 

933 

833 

733 

433 

315 

915 

915 

915 

815 

1 

380ft 

12225  209 

227 

427 

527 

427 

3?7 

227 

109 

309 

509 

509 

409 

1 

3747 

122287501300 

700 

300 

427 

827 

1327 

1327 

1127 

627 

100 

1400 

1500 

1 

3748 

122275000209 

909 

909 

609 

309 

327 

1627 

2327 

2527 

2327 

1627 

727 

1 

3748 

122297501024 

1506 

2806 

3306 

3006 

1906 

124 

1624 

3724 

4524 

4124 

2724 

1 

3749 

122275000508 

2209 

3008 

3108 

2508 

1508 

126 

1026 

1526 

1426 

1126 

526 

1 

3749 

122287500509 

2209 

3009 

3109 

2509 

1409 

124 

1024 

1524 

1424 

1124 

524 

1 

3750 

122275000306 

1506 

2206 

2606 

2206 

1306 

124 

1024 

1624 

1624 

1324 

1024 

1 

3751 

122275000135 

435 

435 

495 

235 

ns 

119 

419 

519 

519 

319 

119 

1 

3761 

122287500100 

100 

100 

100 

100 

100 

109 

109 

109 

109 

109 

109 

1 

375? 

122275000033 

33 

33 

33 

33 

33 

118 

218 

218 

218 

218 

118 

1 

3754 

nr o 

(TOO 

rim 

sim 

fro* 

?m 

m* 

t*tr* 

rw 

rro* 

t 

375ft 

122275001218 

100 

800 

1200 

1200 

1000 

500 

218 

1218 

1818 

2118 

1618 

1 

3  r  sh 

1222/500U318 

104 

304 

4  0  4 

404 

304 

204 

1  18 

318 

518 

618 

518 

1 

3«oo 

1 222750004  1 5 

5?  7 

827 

1027 

1127 

1027 

627 

315 

1015 

1315 

1315 

1015 

1 

3802 

122275000415 

527 

827 

1027 

1127 

1027 

527 

315 

1015 

1315 

1115 

1015 

1 

3804 

122275000119 

200 

400 

bfl  0 

400 

300 

200 

11« 

418 

419 

418 

318 

1 

380* 

122275000409 

27 

27 

27 

27 

27 

27 

209 

709 

1009 

1109 

609 

1 

3744 

12230  400 

90  0 

1000 

900 

600 

100 

4  1  8 

1118 

1318 

1218 

818 

318 

1 

3746 

12230  400 

700 

1000 

900 

600 

100 

422 

1122 

1322 

1222 

822 

322 

1 

3748 

12230  922 

1404 

2604 

3104 

2804 

1704 

122 

1522 

3522 

4322 

3922 

2522 

1 

3749 

122312500922 

1404 

2604 

3104 

2804 

1704 

122 

1522 

352? 

43  22 

3922 

2522 

1 

3749 

12230  709 

209 

627 

1127 

1127 

1027 

827 

827 

827 

209 

909 

1009 

1 

3749 

122312500712 

1612 

1812 

1912 

1412 

812 

112 

727 

927 

727 

527 

127 

1 

3744 

122325000423 

905 

1605 

1505 

1108 

205 

423 

1323 

2023 

2023 

1723 

1123 

1 

3746 

122325000423 

1005 

1605 

1505 

1105 

205 

423 

1323 

2023 

2023 

1723 

1  123 

1 

3749 

122325000427 

1409 

2209 

1709 

1309 

709 

109 

827 

1027 

727 

327 

127 

1 

3744 

12235  121 

903 

903 

1003 

803 

603 

214 

1021 

1721 

1521 

1121 

621 

1 

3746 

12235  121 

903 

903 

1003 

803 

603 

212 

1021 

1721 

1521 

1121 

721 

) 

3749 

12235  121 

903 

903 

1003 

803 

503 

212 

1021 

1721 

1521 

1121 

721 

1 

999 

**r>Er<C.,l) 

8*riEF (m,w,c^OU-TAm,SANFK,01 . CG4 S . DP 0 2 , . ,,«• 1.999999) 
S*OEF  («,  ,CGOIJ-TAM.SANFR,ni  , CGAS . DP02 . UNU5ED > 

f  i 
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LISTING  OF  LITSF 


49 


i 


$J0R,GC0MP62?.  19918FSE.3.  1000 

?SCHFD,C0RE=4n.SC9=9,TIMF=3,CLASS=C.841=l  (16) 
t*OEF  (R.«i,CGOU-TAM.l  ITSF.Ol  .CGAS.0P0  2,  ALL) 

**OEF  (A,  ,CGOU-Tam.LITSF,01  .C5AS.D902.  128  0.400*  ,  .,841  ,  16) 

S*DFF (O.m.LITS.CGOU-TAM.lITSF.OI .CGAS.O) 

*MAP=N 

SFTNU(X=LITS*L> 

PROGRAM  TIOALCUR 

DIMENSION  A ( 37) ♦ AMP (37) .EPOC (37) ,  XODE (114) ,  VPU( 1 14) .MO < 1 3) .SP < 37 > 
1NRDAY (13) » NED A  Y  ( 1  3 ) .XCOS (1025) .SPO ( 37) . AWG<  37) . TABHR (24) ,AKU(S5) ♦ 
PANG (37) ,KDAY (32) .ST09X ( 8  1 * ) , EXT  I M ( 780 ) ♦ JXT I m ( ?60 ) .VEL<260) ,S(37) . 
3EPOCH  (  37  )  .  AMP  A  (37)  .JJXTI<16).JSTIM<260). 

4 I YR (IS) .NUM(IS) *  I  ST  A ( 6 ) , NO ( 6 ) «JJSTI (10) »XVEL(10> 


LTSF 
♦  L  TSF 
LTSF 
LTSF 
LTSF 
LTSF 


COMMON  YVEL(IO) . 

JJJST (10) 

LTSF 

7 

DATA  TABHR/ 

- 

24.  , 

720.,  1392..  2136.. 

2856.,  3600. .LTSF 

8 

1  4320.,  5064. 

«  Sti08«9  65?8*  • 

7272..  7992.,-  24., 

720.,  1416. .LTSF 

9 

2  2160.,  2880 . 

•  3624, 9  4344 • « 

5088.,  5832.,  6552., 

7296.,  8016. /LTSF 

10 

DATA  A«U/ 

7.*, 

7.6, 

7.6.  7.6, 

7.7.  7.7, 

7 

.8,  7.8.  7. ft. LTSF 

11 

1  7.9.  7.9,  8.0, 

8.0,  H.n, 

8.1. 

8.1,  8.1. 

8.2.  8.2, 

8 

.2.  8.3.  8.3.LTSE 

1? 

2  ft. 4,  8.4,  8.4, 

8.5,  8.5, 

8.6, 

8.6,  8,6, 

8.7.  8.7. 

8 

.7,  8.8,  8. 8, LTSF 

13 

3  8.8.  8.9.  8.9, 

9.0«  9.0, 

9.0. 

9.1.  9.1, 

9.2.  9.2, 

9 

.2,  9.3,  9. 3. LTSF 

14 

4  9.3.  9.4.  9.4, 

9.4.  9.5. 

9.5. 

9.6/ 

LTSF 

15 

DATA  */ 

28 

.9841042, 

30.0000000 

• 

28. 4397295. LTSF 

16 

1  15.0410686, 

57.9682084 

.  13 

.9430356, 

86.9523127 

9 

44. 0251729. LTSF 

17 

?  60.0000000. 

57.4238337 

,  ?» 

.5125831 , 

90. 0000000 

9 

27.9682084,1 TSF 

18 

327.8953546.16.1391017,29.4556253,15.0000000.14.4966939, 


4  1 5 .5854433.  0.5443747, 

5  1 .0950331.  13.4715144, 

6  12.8542862.  14.9589314. 

7  42.9271398.  30.0821373,  l 

NNO AY  =  0 

DEVELOP  COSINE  TARLE 
H*. 00153398 

«=2.0-p*H 

MA»T=64 
NAP T=n 

no  35  1=1.16 
NAPT1 =NAPT* 1 
NART2=NART*2 

part=nawt 

phi A=paRT*m 

Pm  I m=pm I  A ,H 

XCOS (NAKT1 ) —COS ( Pm I  a ) 

XCOS (WANT?) »COS (PHI H) 
MAPT1=MAPT-1 
DO  30  ,J=NAWT2.maPT1 
XCOS (J*l)=H*XCOS(J) -XCOS (J-l ) 
30  CONTINUF 
NART=mapT 
maRT=maPT*64 
35  CONTINUE 

XCOS(10?5)=0.f) 

MS*  1 
MY=  1 
MO—  1 

COM=10?4./90. 

DO  90  J= 1 , 37 
A (J) =A ( J) *C0N 
90  CONTINUE 
m(.0  =  0 

95  NSE  0  =  1 

IF  (MS)  120.120.110 
1 1 0  READ  550 


0.0821373, 
13.3986609. 
31 .0158958, 
115.9364169. 


0.0410686. 

29.9589333, 

43.4761563, 

58.9841042/ 


LTSF  19 
1 .01589S8.LTSF  20 

30. 0410667. LTSF  21 
29. 52847 89, LTSF  2? 

LTSF  23 
LTSF  24 
LTSF  25 
LTSF  26 
LTSF  27 
LTSF  28 
LTSF  2Q 
LTSF  30 
LTSF  31 
LTSF  3? 
LTSF  33 
LTSF  34 
LTSF  35 
LTSF  36 
LTSF  97 
LTSF  38 
LTSF  39 
LTSF  40 
LTSF  41 
LTSF  4? 
LTSF  43 
LTSF  44 
LTSF  45 
LTSF  46 
LTSF  47 
LTSF  4ft 
LTSF  49 
LTSF  50 
LTSF  51 
LTSF  5? 
LTSF  53 
LTSF  54 
LTSF  55 
LTSF  56 


Ht  OlJ  1J/  ,  Kt  KMC  •  1  -M U  1  •  1  •  i  qiJ  •»  .  1  ^*.14  ,  1  NIJ5  «  1  NUb  *  I'd-  U  i  K  ,1MtlJ  i  K 

L  1  SF 

5  t 

RE  AO  sn,  ISTA  (  1  )  .N0(  1  >  *  (  A’«P(  J)  ,  F  ROC  (  J  )  .  J=1  •  7  )  ,  1  ST  A  (?)  .NO  (?)  , 

LTSF 

SB 

1  (AMP(J)  ,EBOC  ( J)  .J=R.]4>  .  ISTA(3>  .NO  (3)  ,  (AMP(.J)  ,frOC<J>  .J=15 

.21 )  . 

LTSF 

SO 

?IST<U4>  .NO(*)  .  (AMP(J)  .t-POr(J)  ,J=??.?R)  , ISTA(S)  .NO (6)  *  (AMP(J)  . 

LTSF 

60 

3EPOCIJ)  .  J=?9,  35)  .  1ST  A  <s)  .NO  <b)  .  <AMP(j)  .EPOC(J)  ,J  =  3<S  *  3  7  > 

LTSF 

61 

00  115  L  =  1  » 5 

LTSF 

6? 

IF  (ISTA(L)  .ME.  ISTA(L*D)  GO  TO  451 

LTSF 

63 

115  CONTINUE 

LTSF 

64 

TSTA1  =  ISTA  <1  ) 

LTSF 

65 

00  1 1 R  L  =  1  .  b 

LTSF 

66 

IF  ( NO ( L )  .NE .  L)  GO  TO  450 

LTSF 

67 

1  1  B  CONTINUE 

LTSF 

SB 

120  IF  (MY)  131.131.125 

lt^f 

69 

125  READ  533.  I  Y«  ( 1 )  ,LY  1  *nijm  (  1  )  ,  (  XODE  (  J)  •  VPU  (  J)  .  J=1 .5)  .  I  YR  (2)  . 

LY2. 

LTSF 

70 

1NUM  (2)  .  (  XODF.  ( vJ )  .  VPU  ( J)  » J=9. 16) . IY« (3)  .LY3.NUM  (  3)  .  (XODE  < J)  . 

VPU  (  J )  , 

LTSF 

71 

2J= 1 7  »  ?4 )  . I YR (4)  .LY4.NUM (4)  ,  ( X  ODE ( J) . VPU ( J)  , J  =  ?5 . 32 )  *  I YW ( 5 ) 

.LY5, 

LTSF 

7? 

3NUV<5> . (XODE(J) ,VPU(J) ,J=33.40).IYR<6) «  L  Yb.NlJH ( 6 ) . (XODE(J) 

. VPU ( J) 

.LTSF 

73 

4J=41«4R)  •  IYP  (7)  *LY7«NiiM(7)  ♦  (  XODE  (  J )  , VPU(J)  ,0  =  49.56)  » I YR ( B ) 

•  L  YB  , 

LTSF 

74 

5NUM15)  ,  (XOOE(J)  ,VPU(J)  ,  J=57 .6**)  .  IYP(9)  ,LY9.NUM(9>  ,  (XODE(J) 

,VPiJ(  J) 

.LTSF 

7s 

6 J  =  65  »  72 )  , IYR(IO)  .LYIO.NU'MIO)  .  ( XOOF  (J)  .  VPU(J)  , U= 73  .  BO )  ,  I  YR (  1  1 )  . 

LTSF 

76 

7LY 1 1 .NUM (11). (XOOE ( J)  , VPU ( J)  . J=B1  ,  «3>  . I YR (  1?)  »  L Y 1 2  »  NUM (12) 

* 

LTSF 

77 

B (XODE ( J) .VPU ( J)  . J  =R9.96) « I YR ( 1 3) , L Y 1 3 , NUM (13).  (XODE  ( J) . VPU ( J ) . 

LTSF 

7  B 

9J=97 .104) . I YR ( 1 4 ) ,LY14,NUM(14) . (XOOE(J) .VPU(J) .J= 105. 112) , 

1  YR  ( 1 5 ) 

.LTSF 

79 

1LY15.NHM ( 16) . (XODE ( J) , VPl  MJ) ,J=1 13.114) 

LTSF 

nn 

DO  127  L=1 ♦ 1 4 

LTcF 

R 1 

IF  (IYR(L)  .NF.  IY«(L>1))  GO  TD  45? 

LTSF 

BP 

127  CONTINUE 

LTSF 

BP 

DO  130  L  =  1 • 1 5 

LT-F 

Ba 

IF  <NIJM(L)  .NF.  L>  GO  TO  453 

ITS*-- 

RS 

130  CONTINUE 

l  T  -F 

Bs 

IYP0=moD(IYR(1 ) ,100) 

LTSF 

R  7 

I YR 1 = I YR ( 1 ) 

LTSF 

h  R 

131  IF  (MO)  150.160, 140 

l  TSF 

aq 

140  READ  51*.  (MO ( J) .NBDAY ( J) .NEDAY ( J) , J=1 . 12) 

LTSF 

jn 

SET  UP  TABLES  FOw  N0-7ER0  CONSTITUENTS 

LTSF 

91 

MH0Y=NPOAY ( 1 ) 

M01=M0< 1 ) 

NEDY sNFDA y  (  1  ) 

160  K=0 

LTSF 

OP 

00  1  Bo  1=1,37 

LTSF 

03 

IF  (AMp(J))  1BO.1BO.170 

LTSF 

94 

170  K=K ♦ 1 

lt^f 

*ys 

AMP ft ( K ) =Amp ( j) 4 XODE ( J) 

LTSF 

96 

TEMX=VPU(J)-tPOC( J) 

LTSF 

97 

IF  (TFmx  ,GE.  0.)  GO  TO  171 

LTSF 

9» 

TEMX  =  TEMX  +  3*-0. 

LTSF 

99 

171  EPOCH(K)=TFmx*CON 

LT ‘  F 

mo 

SPD (K)sA(j) 

LTSF 

11 

SP(K)=SPD(<)/10. 

I.TSF 

102 

S(K)=SPO(K)/60. 

ltsf- 

103 

1 b o  CONTINUE 

ltsf 

104 

NOCON=K 

LTSF 

105 

OPERATING  TahlES  NOW  sTOhEO  as  AMP  a  (K  )  .EPOCH  ( *  )  .  SPD  <K  ) 

LTSF 

1  06 

DO  4000  JP= 1 .12 

LTSF 

107 

IF  (MO(JP))  4005.4005, 1S5 

LTSF 

1  0  B 

1«5  M0(13)=MO(JP) 

LTSF 

1  09 

NBDAY (13) aMHOAY ( JP) 

LTSF 

110 

NEOAY ( 1 3) sNEDAY ( JP) 

LTSF 

111 

NNEDA*NE()AY(1  3)  *1 

LTSF 

11? 

NODAYSsNEDAY ( 1 3) -NBDAY ( 1 3) *2 

LTSF 

113 

N0HRS=N0DAYS*?4 

LTSF 

114 

IF  ( MO ( 1 3 )  .NE.  1?)  GO  TO  190 

LTSF 

115 

N0HRS*n0HR5*24 

LTSF 

116 

190  HPS«N0H«S 

LTSF 

117 

OgTgPMfNr  FI  PST  -WJUP  ryf  rTMT  ygPTDO 

xrsr 

1TB 

IF  <LY1)  200.200,210 

LT  bf 

119 

?no  k=mu(1/»> 

GO  TO  PIS 
?l»t  E=MO (131+1? 

PIS  hl»ay  =  mh:)ay  ( )  3) 

PIHST-Tft^H-i  (i*  )  +  80AY*?1*. 

nK  I hST  =f  IHsT 

oo  p?o  jsi.m.s 

STOHX  < J) =0 . 

??0  COM  I  SUE. 

CiJHRKNT=PF»yC  +  rtME>A  (K  )  (  A  (K  )  O  T  *FPOCH  (K  )  ) 

KOUNT  =  n 
K  T  =  0 

00  380  *  =  1  «  N  0  ^  *  S 

If-  (KOUNT  .GT  .  u)  oo  T 0  ?+0 

<<  GUNT  =  1 

?3 1  DO  2 SO  I  = 1 « NOCON 

A8Gu  =  SE,D  i  j  )  *F  I  HST  ♦  EhOC:-4  (  J ) 

AHO  (  J)  =  A  y  0  U  < APGU. 40*6.  ) 

PS n  CONTINUE 
00  TO  ?H0 

?60  DO  280  J=1.N0'~0N 

AhG ( J )  =  AHG  (  J)  +SHf)  (  J) 

P70  IE  (APfi(J)  ,LT.  4096.)  M>  TO  PH) 

A  HO  (  J)  =  A  HG  (  J  )  -409*  . 

GO  TO  ?70 
?8)j  CUMIN.)- 
?90  DO  074  J=  1  .  M)r  ON 

IK  (  ADO  (  J) -1  0?4  .  )  3?<).  J?0.  jl)0 

Ton  Ih  (  AHM  J) -?U4*.  )  Tsn.isn.Tio 
Tin  IE  (  A9G<  J) -T07P.  )  360.380+330 
TPn  AND ( J) =aho ( J) 

00  TO  T40 

T  TO  ANO(J)  =4096. -AhO(.I) 

340  NP=ANG ( J) ♦ 1 .5 

STORX  (K  >  =ST()WX  (IX)  +  AMPA  <  J)  "  XCOS  <NK>) 

00  TO  T  7  4 

TSn  anG(.j)  =204rt.-AHO( J) 

GO  TO  170 

Tftn  ANO ( J) =AHO ( J) -?048. 

T  7n  NHrANO ( J)  ♦ 1  ,S 

STOHX  fK) -STD»X  (K ) -AMP A ( J ) * X COS ( N« ) 

P?4  CONTINUE 

IK  (•<  •  NE  .  NOHHS)  GO  TO  380 
IK  (KT  . K  0 •  1)  GO  TO  3/h 

K  I  HST  =K  I HST  +  HHS- 1 . 

M  =  1 

ChECK-STOhx  <« ) 

STOHX (K ) =0. 

GO  TO  ?  3 1 

T7R  C«SUM=ChECK~STOhx (K) 

380  CONTINUE 

IE  ( I N 0 4  . NK .  1)  GO  TO  398 
DO  38S  K=  1  . NOmRS 
IK  (STOHX(Kl)  3*1.385,38? 

T«1  Nt K  =  0 

STOHX <*) SSTOHX (K) » (-1 .0) 

GO  TO  38 3 
38?  Nfc.F *  1 

3«3  STOHX (*) =SOMT (STOHX (K )  ) 

IK  (NKK  .60.  1)  GO  TO  J8S 
STOHX  (*<■)  =STOHX  O'  )  * (-1 . 9> 

38S  CONTINUE 

395  DO  400  k=i,W0mhS 

STOHX (K ) =STOwX (n ) ♦PEHNC 

4nn  continue 

so  TO  (Alw.4  01 .401 )  .  I'JUl 


LISE 

l?o 

ltse 

1?1 

ltsf 

\?2 

LTSK 

1?3 

LTSK 

1?4 

LTSK 

IPS 

ltsk 

1  26 

ltsk 

127 

ltsk 

128 

ltsk 

IPS 

ltsk 

130 

LTSS 

131 

ltse 

13? 

ltsk 

133 

ltsk 

134 

ltsk 

13S 

ltsk 

1  36 

ltsf 

137 

LTSF 

1  38 

ltsk 

139 

LTSF 

140 

LTSK 

1  4  1 

ltsk 

14? 

LTSK 

143 

ltsk 

144 

LTSK 

14=- 

ltsk 

146 

ltsk 

147 

LTSK 

14« 

LTSK 

149 

LTSK 

ISO 

LTSK 

isi 

ltsk 

is? 

ltsf 

153 

LTSK 

154 

LTSK 

155 

ltsf 

156 

ltsf 

157 

ltsk 

158 

ltsk 

159 

LTSF 

160 

LTSF 

161 

LTSK 

16? 

LTSK 

163 

LTSK 

164 

ltsf 

165 

LTSF 

166 

LTSK 

167 

LTSF 

168 

LTSF 

169 

LTSF 

170 

LTSF 

171 

LTSF 

17? 

LTSF 

173 

LTSK 

174 

LTSF 

175 

LTSF 

176 

LTSF 

177 

LTSF 

178 

LTSF 

179 

LTSF 

1«0 

LTSF 

l«l 

LTSF 

1«? 

LTSF 

183 

LTSF 

184 

LTSF 

185 

401  KIJAT  (  |  )  =  N4L)A  r  (  i  J) 

NODAYS=''iOOAYS-l 
1)0  410  1  =  2*  TODAYS 
KOAY<n=KI)AY<I-l>»l 
410  CONTINUE 
o^INT  440 

PRINT  555,  IY»1,M0(13)  .C«5UM,NF0IR.NEDI* 

PWINT  556 

PRINT  537,  ( KD A  Y  < I )  .STOP*  (  ?4*  I -?  3 )  ,  S  T  OH  x  (  ?4  *  T -2?  )  ,  STOHX  ( 24  *  I -2  I  )  , 
1ST0PX  (?4*I-20)  .STOHa  (2**1- 19)  ,ST3»x (24*1-18) .STOHX (24*1-17) . 
2ST0RX (24*1-16)  .STOP*  <?4*T-15>  .  STOP  x  ( 24*  I  -  )  4  )  .STOHX  (24*1-13)  , 
3ST0PX (24*1-12) «  KOA  Y ( I ) .STOPX ( 24* I - 1 1 )» STOHX (24*1-10)  , 

4ST0PX  <?4*I-9)  ,  STOP  X  (24*I-(i)  .STOHX  (?4*I-7)  .STOP*  (24*1-6)  . 

SSTOPX (24*1 -5 ) .STOHX (24*1-4) .STOP X (24*1-3) .STOHX (24*1-2)  . 

6ST0WX (24*1-1 ) .STOPX (24*1 ) ,  1  =  1 .MOO AYS! 

419  IF  (IN01  .EG.  3)  GO  TO  4000 
ITFMS=0 

EXTIM ( 1 ) =4000. 

K  =  1 
MST=1 

N0HRS=N0HPS-1 
I JOH= 1 

DO  3000  1=1. JOHNS 
GO  TO  (1036.1 055 ) « I JO« 

103s  GO  TO  (1039. 2576.2621. 2641.26 PI ) . NST 
1 0 3P  T IME=NF InsT*1 n 
1040  NrP  = 0 

mwn0A=1 
NAHC=1 
GO  TO  1060 

1050  TIME=(NEIRST*I-2)*10 
MHP  =  0 

10S5  NARC=1 
1 OS0  STOXN=0. 

GO  TO  (1075.1100) .NANC 
1075  DO  1090  I-  1  . •'IOC ON 

IF  (NhP  .FU.  1)  GO  TO  10  7*-, 

A N G U  =  5 P  (  J  )  *  T  I  M F  ♦FMO<  H  ( ,J ) 

GO  TO  liiHP 

1  o 76  ANGU  =  S(J)*TIMF*EP0Ch(.J) 

1  0  8Q  APG  (  J)  =A**OD  (  AHUU.4096.  ) 

1090  CONTINUE 

GO  TO  1 1 2 A 

1100  00  11)0  J=1.YOCON 

IF  ( N,JP  .c-').  1)  GO  TO  11-71 
AHfi(J)  =  AW',.  (  J)  *SH(J) 

GO  TO  1105 

1101  AHG ( J) rAWG ( J) *S ( J) 

1105  IF  (AOG(vJ)  .LT.  4096.)  GO  TO  1110 
APG  (  J)  =AWO  (  J)  -4IJP*-  . 

GO  TO  1105 
1110  CONTINUE 
1120  DO  1220  J= 1 . NOCON 

IF  ( APG ( J ) - 1 024 .  )  1  160.1  160.11  30 
1130  IF  (APG(J)-204«.)  1150.1120.1140 

1140  IF  ( A»G ( J ) -3072. )  1190,1190.1160 

1150  ANG ( J) sAWG ( J) 

GO  TO  1170 

1160  ANG  ( J)  =4096. -A^i,  (  j) 

1170  NP* ANG ( J) ♦ 1  .5 

ST0XW=ST0x9*AMPA (J) *  XCUS (NP) 

GO  TO  1220 

1180  ANG ( J) *204«.-ANG ( J) 

GO  TO  1200 

»N0tU)*A*frTTtt-7rt)*H. 

1200  NP* ANG ( J ) ♦  1  » 6 


LIbt-  1H0 
LTSF  167 
LTSF  14A 
LTSF  169 
LTSF  190 
L  T  4F  191 
LTSF  192 
LTSF  193 
LTSF  194 
LTSF  194 
LTSF  196 
l.TSF  197 
LTSF  ) 94 
LTSF  199 
LTSF  200 
LTSf  201 
LTSF  202 

LTSF  <r  03 
LTSF  204 

LTSS  205 
L  T  Sr  206 
L  T  si-  (-fil- 
LTSF  2A« 
LTSF  2 04 
LTSF  210 
LTsf  2)1 
LTSF  2)2 
LTSF  213 
LTSF  214 
LTSF  214 
LTsf  ?16 
I.Tsf  21  ' 

L  T  sf  2 1 s 
LTS~  214 

LTSF  22" 
LTSF  221 
LT  ^F  222 
LTSF  223 
LTSF  224 
LTSF  224 
LTSF  226 
LTSF  227 
LTSF  ??* 
LTSF  224 
LTSF  230 
LTSF  231 
ITS'  23? 

L’SF  233 

LTSF  234 
LTSF  234 
LTSF  236 
LTSt  ?37 
LTSF  23*> 
L  I  SF  239 
LTSF  240 
LTSF  241 
LTSF  24? 
LTSF  243 
LTSF  444 
LTSF  245 
LTSF  246 
LTSF  247 
LTSF  244 
LTSF  249 
LTSF  250 
LTSF  261 


b  I  Unwsbl  u*K-n *24  l  j)  -XLUb  W2) 
l??0  CONTINUE 

IF  (1*04  .IMF.  1)  GO  TO  1  ??S 
It-  (ST'lXw)  1221,122b,  122  3 
1??1  -JkF  =  o 

STOXP=STOXP* (-1 .0) 

GO  TO  1224 
1 ?P3  NE  F  =  1 

1224  ST0X8=S09T  (STOXH) 

IF  (NFF  .£0.  1)  GO  TO  122b 
STOX9  =  STOXPM-1.0) 

1225  STOX«=STOXP*PF89C 

GO  TO  (2500, 2505, 2510. 2710. 2692, 2684, 2682. 2693, 2715)  .NWhOA 
2*00  nwmOA=2 
'MARC*7 
Pl=STOX« 

GO  TO  1060 
?505  mkH0A  =  3 
P2=STGX« 

GO  TO  10  6  0 
2^10  P3=ST'txP 

IF  (PI)  2515,2520,2525 
2515  IF  (P?)  2530.2535,2575 
2520  IP  (P2)  2b<.5, 2550, 2544 

2525  IF  (“5)  2620,2565.2570 
2530  IF  (P1-P2)  2575.2580.2585 
2539  IF  (P2-P3)  2543,2595.2595 
2543  T  I  m£=T  I  f'E  ♦  1  • 

2^44  jSw=l 

GO  TO  2731 
2545  JSW=2 

GO  TO  2731 

2950  IF  (P3)  2670,2600,26Gb 

2S61  TIvFstlMe*!. 

GO  TO  2S45 

2*65  IF  (P2-P3)  2605, ?60b, 2bbl 
2570  IF  (P1-P2)  26) 0,2625,2620 
?G7S  JS*  =  1 
NSW  =  ? 

NbT  =  ? 

MK)t  =  l 

2^76  IF  (STOWX(I))  2577,2630,2630 
2677  IF  (STORX(I)-STOPX(I»in  2999,2640.2640 
2590  IF  (P2-P3)  2660.2600,2670 
?5»S  IF  (P2-ST0PX  1 1 ♦ 1 > )  267b, 2675,2660 
2S9b  JS#=2 
NSw  =  2 

?S96  TIm£=(TIP£*1.)*6,-3. 

2997  E  X  T I M ( K ) =400  0  . 

FXTIM(K*3) =4000. 

GO  TO  ?90? 

2*00  T  I»»E  =  T  IME* 1 . 

GO  TO  1040 
260b  JS*=1 
NSW  =  1 

GO  TO  2996 

?6]0  IF  ( P2-ST0PX (!♦!))  ?690«?69b,?699 

26?o  JSW  =  2 
f(Sw  =  l 
NST  =  3 
nFOFs? 

2*21  IF  (STOP* ( III  2630.2630,2623 

2623  IF  (STOP*  (  H-STOPX  ( 1*1  )  )  ?700.2Tf)0»2999 

2626  IF  (P?-o3)  2h«S, 2600. 270b 

2630  N*hOAs4 

GO  TO  1050 


LISF  2b2 
LTSF  2S3 
LTSF  254 
LTSF  255 
LTSF  256 
LTSF  257 
LTSF  255 
LTSF  259 
LTSF  260 
LTSF  261 
LTSF  262 
LTbF  263 
LTSF  264 
LTSF  265 
LTSF  266 
LTSF  267 
LTSF  268 
LTSF  269 
LTSF  270 
LTSF  271 
LTSF  272 
LTSF  273 
LTSF  274 
LTSF  275 
LTSF  276 
LTSF  277 
LTSF  276 
LTSF  279 
LTSF  250 
LTSF  ?R1 
LTSF  282 
LTSF  2«3 
LTSF  284 
LTSF  285 
LTSF  286 
LTSF  287 
LTSF  288 
LTSF  289 
LTSF  290 
LTSF  291 
LTSF  292 
LTSF  293 
LTSF  294 
LTSF  295 
LTSF  296 
LTSF  297 
LTSF  ?98 
LTSF  299 
LTSF  300 
LTSF  301 
LTSF  30? 
LTSF  303 
LTSF  304 
LTSF  305 
LTSF  306 
LTSF  307 
LTSF  30R 
LTSF  309 
LTSF  310 
LTSF  311 
LTSF  31? 
LTSF  313 
LTSF  314 
LTSF  31S 
LTSF  316 
LTSF  317 


I 


264H  Nvy-ijfl  =  -i 
JSw  =  2 

2641  tXTIM(K)s4000. 

EXT Im (w+3) =4000. 

2642  IF  ( MO  <  1  3 )  .NF.  I)  6  0  TO  1050 
IF  (I  .ME.  1)  00  TO  1050 

T IME=NF  I PST»  1  n 

mhh  =  o 

GO  TO  10SS 
2660  JSw  =  3 
NSW=1 

2*6 1  TI*F=NFI9ST*60 
GO  TO  29Q? 

2670  JSw  =  2 
NSW  =  2 

TIMF=mF  I9ST»6'' 

GO  TO  2597 

267  =  TIMFrTIMF*^. 

POINT  1 =  °2 
J5w  =  9 
MSW  =  I 

NWM04=6 
I J0B  =  2 
MGw  =  n 

GO  TO  3000 

2680  NWH0A  =  7 
JSW  =  3 
NSW  =  1 
NST  =  4 

2681  IF  <STOWX  (  I  J-STOPX  (  1*1  )  )  ?642. 2642,2999 

?682  NWHOA=S 

26«3  POI NT ] =STOX w 
NAPC=? 

TIMF=TIME+1 . 

GO  TO  1060 

2684  IF  (POINTl -STOXP)  2900.26*6,2740 

2685  JSW=1 
NSW=1 

TIME=NFIPST*60 
GO  TO  2597 
?686  OTIME=TI*F 

TImE=TImE*g.-g. 

GO  TO  2992 

2690  MwMOAsS 
JSW  =  4 

mSW  =  2 
NST  =  S 

269 1  IF  (STOPX ( I ) -ST09X ( !♦! ) )  2999, ?64? , 2642 

2692  NwHOAsA 

GO  TO  2683 

2693  if  (POINTI-STOXK)  2740,2666,2900 
269S  TIME=TIME*2. 

P0INT1=P2 
JSW  =  4 
NSW=? 

N#HOA=» 

IJ0H=2 

NHWsO 

GO  TO  3000 
2700  NwhOA=7 
JSw»l 

GO  TO  2641 
?70S  JSW*4 
NSw  =  2 

"GO  TTT  2G61 
?710  NWH04S9 


LISF  Jl« 
L  TSC  319 
LTSF  320 
LTSF  321 
L  T  SF  322 
LTSF  323 
LTSG  3?4 
LTSF  32S 
LTSF  326 
i.TSF  327 
LTSF  3 pm 
ITSF  329 
LTSF  330 

L  T  ^  f  33, 

LTSF  33? 
LTSF  333 
LTSF  334 
LTSF  33S 
LTSF  336 
LTSF  337 
LTSF  33m 
LTSF  329 
LTSF  340 
LTSF  <4  1 
LTSF  342 
LTSF  343 
LTSF  344 
LTSF  34S 
LTSF  346 
LTSF  347 
I.TSF  948 

LTSF  949 
LTSF  350 
LTSF  3S1 
LTSF  352 
LTSF  353 
LTSF  354 
LTSF  355 
LTSF  356 
LTSF  357 
LTSF  35« 
LTSF  359 
LTSF  360 
LTSF  3hi 
LTSF  362 
ITSF  363 
LTSF  364 
LTSF  965 
LTSF  366 
LTSF  367 
LTSF  363 
LTSF  369 
LTSF  370 
LTSF  971 
LTSF  372 
LTSF  373 
LTSF  374 
LTSF  375 
LTSF  376 
LTSF  377 
LTSF  378 
LTSF  370 
LTSF  380 
LTSF  381 
ITSF  GG? 
LTSF  303 


P01NI  1=510X8 
N)AWC=‘-' 

TiMp=riMr*i . 

GO  To  li^O 

271b  IF  (STGXrt  ,M).  0.)  GO  TO  ?THO 
GO  T0  <2720. ?7?S) .MFOf 
?720  IF  ( S  T  G  X  Ff  .GT.  0.)  GO  TO  27SS 

272]  T!MF  =  TlMf  +  l.  1 
POINT  1  =STOXP 
GO  TO  1060 

2725  IF  (STOXb  .LT.  0.)  GO  TO  ?7SS 
GO  TO  ?7?1 

2710  IF  (Nbp  »EQ.  \)  GO  TO  271b 

273)  E'XTIM(K)sTIPt*fe» 

GO  TO  299b 

271b  F.XTIM(K)=TIMt 
GO  TO  2995 
2740  T IME  =  T IME» 1 • 

POIMT1=GTOXP 
GO  TO  1060 

27bb  IF  <N«P  .f:0.  n)  GO  TO  27*0 
N  6  P  =  0 

GO  TO  273b 
?7b0  NhP=l 

NwhOa=4 

T  I WF  s  (  T  I  M*f  - 1  •  )  *6  * 

GO  TO  lObb 
?O00  OT I ^E=T I MF 

TI’^E  =  TI'J|E*<'.-1. 

2902  AmsO. 

no  29«0  L=  1  •  7 
STOXsT, 

IF  <L  .GT.  1)  GO  TO  201b 

00  2010  Jal.lMOCON 

AMGU  =  b ( J) *T IMF ♦EPOCH  < J) 

APG ( J )  sANOO  <  AS>GU«4  096.  ) 

?9]0  CONTINUE 

GO  TO  2930 

?«lb  00  292b  Jsl.NOCON 
AHG < J) sAPG < J) *S(J) 

?Q?A  IF  (A«G(J)  .LT.  4006.)  GO  TO  ?P2b 

AF-G(J)=AHG(J)-40Ob. 

GO  TO  ?0?o 
2b?b  CONTINUE 
2930  00  ?9b0  Jsl.NOCON 

IF  <A«G(J)  -  1024.)  2 935,29 lb ,2932 
?°32  IF  (AWG(J)  -  2048.)  2443.2943,2931 
?411  IF  <A»G<U)  -  1072.)  2944,2944,2940 
?oiS  ANG( J)sA«G(J) 

GO  TO  ?941 

2940  ANG ( J) =4096 . -APG ( J) 

2941  NPs ANG ( J ) ♦ 1 . S 
STOXsSTOX  +  ANpft  < j) *XCOS (NP) 

GO  TO  ?9S0 

2941  ANG<J>=2048.-AWG(J> 

GO  TO  ?94b 

2944  ANG(J)*A40(J)-204H. 

294b  NP=ANG ( J) ♦ 1 .5 

STOXsSTOX-AMPA  < J)*XCOS (NP) 

?9b<)  CONTINUE 

IF  <INf)4  .NF.  1)  GO  TO  29$b 
IF  (STOX)  29-,)  ,2955,2963 
?9b 1  NEF*0 

ST0X*ST0X«<-1.) 

29*?  3  NF. F *  1 

GO  TO  2454 


LIbF  JH4 
LT5F  3bb 
LTSF  386 
LTSF  387 
LTSF  388 
LTSF  189 
LTSF  390 
LTSF  391 
LTSF  34? 
LTSF  391 
LTSF  394 
LTSF  39b 
LtSF  396 
LTSF  397 
LTSF  398 
LTSF  ■•99 
LTSF  400 
LTSF  401 
LTSF  40? 
LTSF  403 
LTSF  404 
LTSF  40b 
LTSF  406 
LTSF  407 
LTSF  4^8 
LTSF  409 
LTSF  4  1  o 
LTSF  411 
LTSF  -1? 
LTSF  413 
LTSF  414 
LTSF  4)5 
LTSF  416 
LTSF  417 
LTSF  418 
LTSF  419 
LTSF  420 
LTSF  421 
LTSF  422 
LTSF  423 
LTSF  424 
LTSF  42b 
LTSF  426 
LTSF  427 
LTSF  428 
LTSF  429 
LTSF  430 
LTSF  431 
LTSF  43? 
LTSF  411 
LTSF  414 
LTSF  43b 
LTSF  436 
LTSF  417 
LTSF  436 
LTSF  439 
LTSF  440 
LTSF  441 
LTSF  442 
LTSF  441 
LTSF  444 
LTSF  445 
LTSF  446 
LTSF  447 
LTSF  449 
LTSF  448 


?9S4  5I0X  =  S<JK|  tSIUXI 

IF  (NFF  .£0.  ])  go  TO  ?9ss 
ST0X=ST0XM-1.) 

?955  STOX=STOX*PE«mC 

IF  (L  .EO.  1)  SAVIT=STOX 
GO  TO  (2960.295b) ,NS* 

2956  IF  (SAVIT  .GE.  STOX)  GO  TO  2980 
GO  TO  2 962 

29*0  IF  (SAVIT  .LE.  STOX)  GO  TO  2940 
2962  SAVIT=STOX 
AM=L-1 

2950  CONTINUE 

EXTIM(<*1)=TIMF*AM 

EXTIM (*♦?) =SA V 1 T 

IF  (K  . r. 0,  1)  GO  TO  2990 

IF  (EXTIM  ( K  ♦  1 )  .GT.  EXTIM  (K-?)  )  GO  TO  2.990 
TIME=0TIME*1  . 

POINTl=STOX» 

GO  TO  10SS 
2990  K  =  K»3 

ITEMS= I TEMS+3 

299S  GO  TO  ( 2690  *  2660 « 25  76* 26?0 ) .  JSX 

2999  IJ06=l 

3000  CONTINUE 

IF  (IMOS)  3081 .3051.3075 
307S  J1TEM=ITEms-3 

DO  30«0  K=1 « JITEM.3 

IF  (EXTIM(K*?)  .GF.  .OS)  GO  TO  30pO 
IF  (FXTIM(K*2)  .LE.  -.05)  GO  TO  3080 
EXTIM (K ) =4000. 

EXTIM(*f*3)=4000. 

30"0  CONTINUE 
30A1  KAY=ITFMS/3 

J  =  0 

DO  3 1 SO  K  =  1  . ITEMS*  3 
r  L  =  « 

J=J*1 

M=f| 

IF  (EXTIM(KL)  .NE.  4000.)  GO  TO  310S 
JSTIM(J) =4000 
GO  TO  3106 

3105  jh8S=rxTlM(i<L) 

JOA  YsMOO (JHPS.  1440) 

JH«  =  Jf'A  Y  /SO 
JMlN=MOD ( JOF Y.60) 

IF  (v  . ) Q.  1)  GO  TO  3110 
JST IM ( J) = JHR*1 00+ JMI N 

3106  KL=KL*1 

M=1 

GO  TO  31  OS 

3110  JXTIM(J)=JmU*ioO*JMIN 
VEL( J)=FXTIM(K»2> 

IF  (VFL(J))  31 1S.3150.31SO 
31 1 S  IF  (VFL(J)  .LF.  -.05)  GO  TO  31SO 
VFL ( J) *VEL ( J) * (-1 . ) 

3150  CONTINUE 
KKseIXAY 

IF  (JP  .EO.  1)  GO  TO  3170 
IF  (JXTIM(I)-NSAV)  31SP.315S.3160 
315?  5A*KAY-1 

00  3153  K0=1  .M 
JSTIM(K0)*JSTIm(K0*1 ) 

JXTIM(KO)=JXTIM(KOM  > 

VkL(K0)*VEL(X04l) 

5195  CONTINUE 

GO  TO  3170 


L  1  SF 

450 

LTSF 

451 

LTSF 

45? 

LTSF 
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LTSF 
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LTSF 

455 

LTSF 
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LTSF 

457 

LTSF 
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LTSF 

459 

LTSF 

460 

ltsf 

46) 

LTSF 

46? 

LTSF 

463 

LTSF 

464 

LTSF 

465 

LTSF 

466 

ltsf 

46  7 

ltsf 

465 

LTSF 

469 

LTSF 

4  7  r, 

LTSF 

47) 

LTSF 

47? 

LTSF 

47? 

LTSF 

474 

lts^ 

475 

LTvf 

476 

LTSe 

477 

LTSF 

475 

I.TsF 

479 

LTSF 

45  0 

LT--.F 

48) 

LTSF 

45? 

LTSF 

45? 

LTSF 

454 

LTSF 

4  P  S 

LTSt 

486 

LT  '-F 

487 

LTSF 

455 

LTSF 

4H» 

LTSF 

490 

LTSF 

401 

LTSF 

49? 

LTSF 

4'J3 

LTSF- 

494 

LTSF 

495 

ltsf 

496 

LTSf 

447 

LT  -  f 

498 

LTSF 

499 

L  T  SF 

son 

I.TSF 

501 

LTSf 

-.0? 

ltsf 

503 

LTSF 

504 

LTSF 

SOS 

LTSF 

S06 

LTSF 

so  7 

LTSF 

608 

LTSF 

5  0  Q 

LTSF 

610 

LTSF 

611 

LTSF 

51? 

LTSF 

513 

wTSF 

51* 

lits 

615 

Jib's  Jb  I  14(1)  *NbAV5 
GO  TO  3170 

3160  J5TIM(KK*1)=JSTIM(KK) 

JXTIM(KK»1 )=JXTIM(KKI 
VEL(K*»1)=VFL(KK) 

KKsKK-1 

IT  (KK  .E o.  0)  GO  TO  3166 
GO  TO  3160 

3165  JST  IM  ( 1  )  3*00  0 
JXTIMd  )=NSAV 

veui>=SAv 

3170  NDAY=N80AY<13) 

NCOUNT=0 
NNJ=1 
k  =  i 

GO  TO  <3174,3174,3178) , IN02 
3174  Phi  NT  550 

PPINT  555, I YH 1 «MO  < 1 3 ) .CKSUM.NFOIS.NEDI* 
PPTNT  675 

31 7H  IF  ( I NOG  .EO.  1)  GO  TO  3165 
<AYK=KAY-1 
OO  3190  Jsl.KAYK 
IF  (VFl.(J)  .Gy.  .25)  GO  TO  3160 
IF  ( VCL  < J)  .LT.  -.2b)  GO  TO  31«0 
JST  IM ( J ) =4000 
JSTIM(J*1)=4G00 
31 «0  CONTINUE 

3186  IF  (IM06)  3200.3200.3166 

3166  DO  3193  1=1. kay 

IF  ( VCL  ( I )  .GF.  12.95)  GO  TO  3 1  •=*  1 
IF  (VSLd)  .LE.  -12.95)  GO  TO  3192 
IF  (VFL(I)  .GE.  7.55)  GO  TO  3167 
IF  ( VrL  < I )  .GT.  -7.55)  GO  TO  3193 
TEPPsVEL ( I ) * (-1 . ) 

LPN  =  1 

nO  TO  3166 
31 6 y  TEMPsVELII) 

LPN=2 

si«6  N=TEMO*10.0-74.5 

GO  TO  ( 3169,3190) .LPN 
3169  VfcL<I)*AKU(N)*<-l.) 

GO  TO  3193 

3190  VELd)=AKU(N> 

GO  TO  3193 

3191  VEL ( I ) =9.6 
GO  TO  3193 

319?  \/EL(I)=-9.6 

3193  CONTINUE 

3200  DO  3350  1=1. KAY 

IF  (JSTIN(I)  .LT.  4000)  GO  TO  3215 
3210  IF  (JXTIMd)  .GT.  JXTIM(Idl)  GO  TO  3250 
JJST I (K ) *9999 
JJXTI  (K)aJXTlMm 
XVEL (* ) =YEL ( I ) 

GO  TO  3310 

3215  IF  (JXTIM(I)  .LT.  JSTIMd))  GO  TO  3210 

IF  (JXTIMd)  .GT.  JXTIM (1*1)1  GO  TO  3230 
JJSTI (K)=JST1M(1> 

JJXTI <K) =JXT IM( II 
XVEL<K)=YEL(I) 

GO  TO  3310 

3219  JJSTI IK) *JST1« ( I ) 

JJXTI (K 1 bJXT Im ( 1 ) 

XYEL  IK ) *VEL ( I ) 

IF  ( JST IM ( I ♦ 1 )  «EQ.  4000)  GO  TO  3260 
IF  (JXTIMd)  .GT.  JSTIM(T*1))  GO  TO  3260 


LI  lb  519 
LTTS  520 
LTTS  521 
LTTS  522 
UTS  523 
LITS  524 
LTTS  525 
LTTS  526 
LITS  527 
LTTS  526 
LITS  529 
LTTS  530 
L'TS  531 
LITS  532 
LTTS  533 
LITS  534 
LTTS  535 
LTTS  536 
LTTS  537 
LITS  538 
LTTS  539 
LTTS  540 
LTTS  541 
LITS  542 
LITS  543 
LTTS  544 
LITS  546 
LTTS  546 
LITS  547 
LTTS  546 
LITS  549 
LITS  550 
LITS  551 
LITS  55? 
LITS  553 
LITS  554 
LITS  565 
LITS  556 
LITS  557 
LITS  556 
LITS  559 
LITS  560 
LITS  561 
LITS  56? 
LITS  563 
LITS  564 
LJTS  565 
LITS  566 
LITS  667 
LITS  568 
LITS  569 
LTTS  570 
LITS  571 
LITS  57? 
LITS  573 
LITS  574 
LTTS  575 
LITS  576 
LITS  577 
LITS  578 
LTTS  579 
LITS  580 
LTTS  581 
LITS  58? 
LITS  583 
LITS  584 


UU  ID  J.100 

LI  IS 

586 

3?sn 

IF  ( JST I M ( I ♦ 1 )  .EQ.  4000)  60 

TO  3?6S 

LITS 

68b 

IF  (JXTIM(I)  ,6T.  JST I M ( I ♦ 1 ) 

)  GO  TO  3290 

lits 

687 

GO  TO  3296 

LITS 

688 

3255 

JJSTI <«)=99 99 

LITS 

689 

JJXTI (K)sJXTlM(I) 

LITS 

690 

XVEL(<)=VEL(I) 

LITS 

691 

3260 

nlast=mnj*ncount 

LITS 

69? 

GO  TO  (326*. 3?b*. 3265) . IND2 

LITS 

593 

3?64 

PRINT  585.  NUAY. (JJSTI (J> .JJXTI (J> .XVEL(J) .J=l.K) 

LITS 

694 

IF  (IFL6  .EQ.  1)  GO  TO  7602 

LITS 

596 

I  t<c= i 

lits 

696 

760? 

00  7600  U=1.6 

LITS 

597 

JJJ6T ( I 6C ) = JJST I (L) 

lits 

698 

YVEUIHC)sXVEL(L) 

LITS 

599 

IBC=I6C*1 

LITS 

600 

7600 

CONTINUE 

LITS 

bOl 

7604 

iflg=i 

lits 

bO? 

IF  (iNn?  .EQ.  1)  GO  TO  3276 

LITS 

603 

3266 

JK=K*1 

LITS 

604 

DO  3266  N=JK,6 

LITS 

606 

JJSTI (Nl=9999 

LITS 

606 

JJXTI (N) =9999 

LITS 

607 

XVEL(N>=99.9 

lits 

bOM 

3  266 

CONTINUE 

LITS 

609 

IF  ( I NO 3  .ECO.  3)  GO  TO  3267 

LITS 

bln 

IF  ( I ND3  .EQ.  1)  GO  TO  3266 

LITS 

611 

3267 

CONTINUE 

LITS 

612 

3268 

NSEQ=MSEQ*1 

LITS 

613 

IF  <K  .GT.  6  .AND .  1ND2  .EQ. 

3)  GO  TO  3288 

lits 

6)4 

3?75 

NNJ=NLAST ♦ 1 

LITS 

615 

NCOUNT=  0 

lits 

bib 

K  *  1 

LITS 

b  1  7 

NOAY=NDAY*l 

LTTS 

bl  8 

NNOAY=NND A Y ♦ 1 

LITS 

6)9 

IF  ( NO A  Y  .NE.  NNEOA)  GO  TO  3360 

LITS 

620 

IF  (MO ( 13)  .Nc .  12)  GO  TO  3287 

LITS 

621 

IF  (NO*Y  .NE.  3?)  GO  TO  3287 

LITS 

62? 

PRINT  S60 

LITS 

623 

PPINT  665.  IYP1.ISTA1 

LTTS 

624 

PPINT  680 

LITS 

626 

PPINT  686.  NDAY. (JSTIM(J) ,JXTIM(J) ,VEL(J) .J=NNJ.KAY) 

LTTS 

626 

3287 

NSAVS=JSTIM(IM  ) 

LTTS 

b?7 

NS A V® JX  T I m ( I ♦ 1 ) 

LITS 

628 

SAVsVFL (1*1 ) 

LITS 

6?9 

GO  TO  4000 

LITS 

630 

3288 

PPINT  570 

LITS 

631 

NNLASsNLAST*! 

LTTS 

63? 

PRINT  596,ISTA1.M0(13) «N0AY. 

IYR1. (JSTIM(J) ,JXTIM(J) .VEL(J) . 

LITS 

633 

1  JsNNJ.NNLAS) 

LITS 

634 

GO  TO  3275 

LITS 

636 

3290 

JJSTI (K ) *9999 

LITS 

636 

JJXTI (K)xjxTIm(I) 

LITS 

637 

XVEL (K) =  VEL ( I ) 

LITS 

638 

GO  TO  3260 

LITS 

639 

3295 

JJSTI (K)s9999 

LTTS 

640 

JJXTI (K)*JXTIM(I) 

LITS 

6*1 

XVEL(«)=VEL(I> 

LITS 

6*P 

3300 

5*K  ♦  1 

LITS 

6*3 

JJST I (K ) *JSTIm ( I ♦! ) 

LITS 

6** 

JJXTI (K ) *9999 

LITS 

6*6 

XVEL (K) s99,9 

LITS 

6*6 

GO  TO  3260 

LITS 

6*7 

3310 

NCOUNT*NCOUNT ♦ 1 

LITS 

6*8 

CTTS 

6*9 

3360 

CONTINUE 

LITS 

650 

3350  CONTINUE 


40(10  CUN  I  1 'Hit  LI  lb 

40(m  MEAD  SI*, MS, MY, MO  L  T  T  S 

IF  (MS*MY*MP>  402D«40?0.os  LTTS 

4 ci?n  CALL  SLACK (JJJST.YVEL, NNOaY.NHOY. MO, IY-d  ,NEPY)  LITS 

IFINNOAY  ,NF ,  0)  GO  TO  4020  LITS 

STOP  LTTS 

450  PR I  NT  SOI  LITS 

STOP  LITS 

4 5 1  PRINT  SO?  LITS 

STOP  LITS 

4S?  PRINT  SO 3  LITS 

STOP  LITS 

453  PRINT  S 0 4  LITS 

STOP  LITS 

sni  F0RhaT(27h  STATION  CAPOS  OUT  OF  ORDER)  LITS 

SO?  F0PMATI31H  STATION  NUMriEt>S  NOT  CONSISTENT)  LTTS 

503  F0R4ATI2MH  YEAR  NUMhEPS  NOT  CONSISTENT)  LITS 

504  FOPMA T ( 24H  Y£ak  CAPUS  OUT  OF  ORDER)  LITS 

S3)  FORMAT (?I4,FS.3.F4. 1,F5.3,F4. l,FS.3tF4. 1 ,F5.3.F4.1 .E5.3.F4.1 ,  LITS 

1FS.3.F4.1,FS.3.F4.1 )  LITS 

S3?  FORMAT (FS.3.6I2.2I4.54X )  LITS 

s  33  FORMAT  < I4.212,f 4.3.F4. 1 ,=4.3«E4. 1 .E4.3.F4. 1 .E4.3.F4. 1 •  LITS 

1F4.3*F4.1,F4,3«F4.1,F4.3.F4.1.F4.3*F4.1)  LITS 

S34  Format ( 36I?*8X)  LITS 

S 3 Y  EO»MAT(IR.l?ER.l)  lits 

538  F OR  1 A T (314)  LITS 

SSO  F ORMAT ( 80H  LITS 

t  LITS 

SSS  FORMAT ( /3?h  PREOICTtO  TIOAL  CUPPENT  YEAR. 15. 8H  MONTH, 13. 1 2H  LITS 
1  CHECKSUM. FI?. 7/17-1  FLOOD  0 1  RFC  T  ION  ,  ,  1 4 , 7H  TRUE . . 1 5X 1 8HEBR  (-)  DLITS 
?I RECTI ON., I4.7H  TRUE./P8H  NO A A ,  NATIONAL  OCEAN  SURVEY/)  LITS 

SSS  FORMAT  (?7h  HOi/kLY  VELOCITIES  IN  KN0TS/7X  111  HD  AY  HOURS  HOURS  LITS 
1  HOURS  HOURS  HOURS  HOURS  HOURS  HOURS  HOURS  LITS 
?POUBS  HOURS  HOUHS/14X103H0/1?  1/13  2/14  3/1S  LTTS 

3  4/lft  S/17  S/18  7/19  8/20  9/21  10/22  11/LITS 

4?3/ )  LITS 

SSS  FORMAT ( /3SH  PREDICTIONS  REGINNING  DECEMBER  32, » IS , 5X 1 OHSTA .  NO.  ,  LITS 

1 1 4 / )  lits 

=  70  FORMAT  I///23H  TROUBLE  DAY  FOLLOWS  /)  LITS 

S7S  format (SX1 1 1HSLACK  MAXIMUM  SLACK  MAXIMUM  SLACK  MAXIMijLITS 

1M  SLACK  MAXIMUM  SLACK  MAXIMUM  SLACK  M AX IMUM/ , SX 1 1 1 HL I TS 
?W ATER  CURRENT  WATER  CURRENT  WATER  CURRENT  WATER  LITS 

3CURRF.NT  wATFR  CURRENT  WATER  CURRENT/1 18H  DAY  TIME  TIMELITS 

4  VELOC  TIME  TIME  VELOC  TIME  TIME  VELOC  TIME  TIME  VELOC  LITS 

STIME  TIME  VELOC  TIME  TIME  VEL0C/118H  H.M.  H.M.  KNOTS  LITS 

SH.M.  H.M.  KNOTS  H.M.  H.M.  KNOTS  H.M.  H.M.  KNOTS  H.M.  H. ML  ITS 

7.  KNOTS  H.M.  H.M.  KNOTS)  LITS 

Sflo  FORMAT  <11 9H  DAY  TIME  TIME  VELOC  TIME  TIME  VELOC  TIME  TIME  LITS 
l  VELOC  TIME  TIME  VELOC  TIME  TIME  VELOC  TIME  TIME  VELOC)  LITS 

505  FORMAT(lH0»I3»6(lG»16»F7.2)/4X.t>(I6»I6,F7.2)/4X,6(I6»IG.F7,2) )  LITS 

800  FORMATd?,  14,312. 5114, 14, FS.l)  .13)  LITS 

SOS  F0PMAT(15H  ST 4 •  NO.  .I4.10H  MONTH  .I2.8H  DAY  ,I?,9H  LITS 

1  YEAR  .I4//BX6(I6.I6.F7.) ) /BXG < 16. IS.E7. I > )  LITS 

END  LITS 

SUBROUTINE  SLACK  ( JSLAK ,CEL .NNOAY . N3DY .MO . I YR 1 ,NEDY )  LITS 

DIMENSION  CUR(I2)  .MOIRU?)  ,JHR(10)  .JMIN(IO)  LTTS 

OIMENSION  JSL AK ( 1 0 ) » CEL  1 1 0 )  LITS 

INTEGER  DAYCONT  LITS 

C  CUR (12)  -  CURRENT  SPEEDS  IN  KNOTS  ANO  LITS 

:  MO I R ( 1 2 )  -  ASSOC  I A  TEO  CURRENT  OIRFCTIOnS  IN  DEGHEES  TRUE  FOR  LITS 

C  1/12TH  OF  TIDAL  CYCLE. 6  EACH  FOR  EBB  ANO  ELOOD.  LITS 

Z  AT  OESIREO  LOCATION  IN  POSITION  MATRIX.  LITS 

Z  JSLAK (10)  -  TIMES  OE  SLACK  WATER  < /ONF  TIME)  FOP  EITHFP  ONE  OR  LITS 

C  TWO  DAYS  OF  TmE  YEAH  (PASSED  From  main  PROGRAM),  LITS 

C  CEL (10)  -  CURRENT  SPEEDS  IN  KNOTS  ASSOCIATED  WITH  TIMES  OE  MAXIMUM  LITS 

C  FLOW.  LITS 


JMM(IO)  AND  JMIW(IU)  -  A-MWAYS  UStO  IN  CALCULATION 

OF  NIJMbEH 

UF 

HOURS 

LI  IS 

717 

OP  OCCURRENCE  8EF08E  Oh  AFTFh  CLOSEST  SLACK 

ft  A  TFh 

T  T  MF 

(JSLAK ) 

LITS 

710 

HEAD  (*0.797)  MONLIST 

LITS 

719 

Rf.AO  (*0.797)  DAYCOwT. IOTSC 

LITS 

720 

797  P 09m AT  (11. 1H, ID 

LITS 

721 

IF  ( D A V CONT . F Cl , 2 )  Nt>Or  =  Nfor 

LITS 

722 

RFAO  POSITION  OF  OCCURRENCE 

LITS 

723 

IF  (IOISC.EO.l )  00  TO  <*19*t 

LITS 

724 

WEAD(*0.13)  ILAD. ILAM, ILOO, ILO^ 

LITS 

7?S 

l\  POPAAT  (I?.  12.  IX,  13.14) 

LITS 

72* 

IF  ( I L AD  . F 0 .  99)  GO  TO  9a 

LITS 

7?7 

GO  TO  4 DM9 

LITS 

720 

4199  REAO  (AO. 401D  ILAD.  ILAM,  ROD,  ILOM 

LITS 

729 

401T  FORMAT ( I?. I?. 13. 12) 

LITS 

730 

IF  (ILAD  ,E0.  99)  GO  TO  94 

LITS 

731 

READ  TIME  AND  DATE  OF  OCCURRENCE 

LITS 

732 

4089  REA 0 ( *0 . 1 )  NhR.NmIN 

LITS 

733 

1  FORMAT < 12. 1?) 

LITS 

734 

IF  (MONLIST. EG. 1 )  GO  TO  93 

LITS 

7  35 

NT  I  ME sMHR* 100  *NM In 

LITS 

73* 

00  8  K=1,10 

LITS 

7  m  7 

JHP (K) = JSL AK ( K ) /100 

L  T  TS 

730 

9  JMIN(K)=JSLAK(K)-(JhR(K)*100) 

LITS 

739 

IE  tDAYCONT.EO.  1  )  NNf)AY=l 

LITS 

74  'i 

IF (0AYC0NT.E0.2)  NML>A Y  =  ? 

LITS 

7  41 

K=G*NMOAY-tt 

L  ITS 

74? 

DETERMINE  TIME  OF  SLACK  kATFR  CLOSEST  TO  TIME  OF 

OCCURRENCE 

LITS 

74? 

5  IF  ( JSLAK <K ) .Eu.9999. A  NO. NT  IME.GT. JSLAK (K 

♦1)  ) 

GO  TO 

800 

LITS 

744 

IF  (NT  I MF  .GT.  JSLAK(K)  .ANO.  NTIMF  .LT. 

JSLAK  ( K  ♦  1  ) 

.ANO. 

JSLAK 

LITS 

745 

1  ( K  ♦  1  )  .NF.  9999)  GO  TO  7.ip 

LIT* 

74* 

IF  ( JSLAK ( K )  ,EO.  99m9  .AND.  NT  I  ME  .LT. 

JSLAK  (K*D)  '-0 

TO 

700 

LITS 

747 

IF  (NT  I  ME  ,  GE  .  JSLAK (K)  .'NO.  JSL*K(K*1) 

.FO. 

9999) 

GO  TO 

700 

LITS 

748 

IF  (NTIMF  .GT.  JSLAK(K)  .AND.  JSLAK(K*1) 

.ED. 

0)  GO 

TO  700 

LITS 

749 

IF  (JSLAK  (K)  .GT. NTI«F  i  >..'•)  TO  700 

UTS 

7  S  n 

GC  TO  *00 

LITS 

7*1 

700  IF  (JSLAK(K)  .(-j. 999». OR.  JSLAK  (K)  .GT. NTIMF)  K  =  K 

-I 

LITS 

75? 

IF  (JSLAK(K) .PO.ni  K=K-1 

LITS 

753 

KKK  =  K 

LITS 

754 

KSL  AK= JSL AK (K) /100 

L  ITS 

755 

XSLAK~KSLAK«10D. 

L  ITS 

7** 

YSLAK= JSLAK (K) -*SLAK 

LITS 

7*7 

ZSLAK=YSLAK/mO. 

LITS 

7*  « 

TSLAK=KSLAK./9LAK 

LITS 

7-9 

IF  (JSLAK (K+i) .to. 9999)  K  =  K*] 

LITS 

74n 

IF  (JSLAK ( K ♦ 1 ) .  t  Q  •  0  )  GO  TO  701 

IDS 

7*1 

KSL  AKr  JSLAK  (  K  ♦  1  )  /  1  0  ft 

LIT* 

7*? 

XSLAKr*SLAK*10f>. 

LTTS 

7*  ? 

YSLAKr JSLAK ( K ♦ 1  I - XSL  AK 

1  TS 

7*4 

7SL AK=YSLAK/60. 

LITS 

7  ft  * 

TTSLAK=KSLAK./SLAK 

LITS 

7  ** 

IF  ( JSLAK (K+l ) .LT. JSLAK (K) )  TTSL AE=TTSL AK*?a . 

LITS 

7*7 

70?  tsran=ttslak-tslak 

LIT* 

7*8 

GO  TO  7i)3 

L  TTS 

7*9 

701  IF  < JSLAK (K*2) ,EO. 9999)  K=k+1 

LITS 

770 

KSL AK= JSLAK (K*2) /100 

LITS 

771 

XSLAK=KSLAK*10O. 

LITS 

77? 

YSLAK= JSLAK ( K ♦? ) - XSL AK 

LITS 

773 

ZSLAK=YSLAK/*o. 

LTTS 

774 

TTSLAKsKSLAK*ZSLAK 

LITS 

775 

TTSLAKsTTSLAK.24, 

LITS 

77* 

GO  TO  702 

LTTS 

777 

703  percol=tspan/*. 

LIT* 

778 

GO  TO  20 

LITS 

779 

"00  K*K ♦ 1 

LITS 

780 

If  fX.TT .  .AfNOAT-*^)  )  8(7  m  70 

-L1TS 

781 

GO  TO  5 

LITS 

782 

20  *  =  L  1  I  S  /HI 

X  M I M:NM I N/ 6  0 .  LITS  784 

XNTIMC-  =  NH(3*>MrN  LITS  785 

IF  (  JSLAK(K)  .GT.NTlMf  )  XNT  I  ME  =  XNT  I ME  ♦  ?(.  .  LTTS  786 

hRaFT=XNTImE-TSLAk  LITS  787 

XJJ1=HRAET  LITS  788 

determine  factor  ry  which  current  spefos  at  desired  location  are  lits  789 


HE 

MULTIPLIED  TU 

GET 

ACTUAL  CURRENT 

L'TS 

790 

IF 

(CEL (<) .Eg. 99 

.9) 

KsK*  1 

LTTS 

7  9  1 

IF 

( CrL  IK)  .Fu.O. 

)  K  = 

K  ♦  1 

LITS 

792 

IF 

(CFL(K)  .GT. 

0.  ) 

GO  TO 

122 

LITS 

793 

IF 

I  IOISC.EO.P) 

GO  TO  600 

LIT0 

794 

SECTIONS  FOR  THE 

PACE  LONG  ISLAND  SOUND 

LITS 

79S 

IF 

<AHS(CEL(K) ) 

.GE. 

1.8 

.  ANO. 

AHS (CEL <*  >  ) 

.LT. 

2. 2) 

FAC=  .5 

LITS 

796 

IF 

( ARS (CEL ( K ) ) 

.  GE  . 

2.2 

.  ANO. 

AHS  (CEL  (K ) ) 

.LT. 

2.6) 

FAC=  .6 

LITS 

797 

IF 

(ARS (CEL (K) ) 

•  GE. 

2.0 

.AND. 

AHS (CEL (K) > 

.LT. 

3.0) 

FAC=  .7 

LITS 

79B 

IF 

(ARS (CEL (K) ) 

.GE. 

7.0 

.  AND. 

AbS (CEL  (« ) ) 

.LT. 

3.4) 

FAC=  .0 

lits 

799 

IF 

(  A  8  S  (  C  E  L  ( K  )  , 

.  G£ . 

3.4 

•  4ND  . 

A8S  (CEL  (<  > ) 

•  LT. 

3.8) 

FAC=  .9 

LITS 

000 

IF 

(ARS (CEL (K) ) 

.GE. 

3.0 

.AND. 

AHS ( CEL  £  K ) ) 

.LT. 

4.2) 

FAC=1.0 

UTS 

ROl 

IF 

(ARS (CEL  (K) ) 

.GE. 

4.2 

.AND. 

ABS  (CEL  (>•  )  ) 

.LT. 

4.6) 

FAC=1 . 1 

LlTv 

«02 

IF 

( ARS I  CEL ( K ) ) 

.GF. 

4.6 

.AND. 

AHS  (CEL (K ) > 

.LT. 

S.O) 

F  AC  =  1 , ? 

LTTS 

003 

IF 

( ARS ( CEL  ( K  )  ) 

•  GE. 

S.O 

.  AND. 

AHS ( CEL ( K ) ) 

.LT. 

S.4) 

FaC=1 . 3 

LITS 

804 

GO 

TO  4 

lits 

80S 

=>00  CONTlNUF 

C  CORRECTIONS  FOR  THE  GOLDEN  GATE 

IF  (ABS(CFL(K)  )  .GT. 0.0. AND.  ABS ( CEL ( < ) )  .LT  .  0 . 7 ) 
IF  (ARSICELCO  )  .GF.0.7.ANO.  A8S  ( CEL  ( K  )  )  .  L  T  .  1 . 2 ) 

IF  (AHSICFL(K)  )  .GE.l  .2. AMD.  AHS  ( CEL  ( K  )  )  .  L  T  .  1  .  6 ) 

IF  (ARS(CEL(K) ) .GE. 1.8. AND.  A8S(CEL(K) ) .  LT . 2 . 1 ) 
IF  (AHSICELfd ) .GF.2.1.AWO.  AHS  (CEL (K ) ) .LT.?.S) 
IF  (APS (CEL (K)) .GE.2.S.AMD.  ARS < CEL < K )) . LT . 3 . 0 > 
IF  < A»S (CEL <  K )  ) .GE. 3.0. AND.  AbS ( CEL ( * ) )  . L T . 3 . 4 ) 
IF  (  APS (CEL C  K  > ) .GE.3.4.AM0.  ARS (CEL (K > ) .LT.3.9) 
IF  <4«S(CEL(K) )  .GF.3.9.AND.  AHS ( CEL ( * ) ) .LT . 4 . 3 > 
IF  < 4«S (CEL (K ) ) .GE.4. 3. AND.  ARS < CEL < K > ) . L T . 4 . 8 ) 
IF  (AHS (CEL (Kl  ) .6E.4.8.ANU.  AHS (CEL ( K )  )  .LT.S.2) 
IF  (AHSICEL (K) ) .GF.S.2.AND.  ARS ( CEL ( < ) ) . LT . S . 7 ) 
IF  (A«S (CEL <K)  ) .GE. 0.7. AND.  AHSICEL (*)  )  .  L  T  ,  6 . 1  ) 

IF  (ARS(CEUK)  )  .Gt". 6.1. AND.  ARS  (CEL  <K)  )  •  L  T  .  6  •  6  ) 

IF  (AHS  (CEL  (r  )  )  .GE.6.6.  AND.  ABS(CELIK) ) .LT.7.1) 
GO  TO  14 

12?  IF  tiOISC.EO.2)  GO  TO  SOI 
C  CORRECTIONS  FOP  THE  RACE  LONG  ISLAND  SOUND 


SOI  CONTINUE 

C  CORRECTIONS  FOR  TrtE  GOLDEN  GATE 

IF  ( ARS  < CEL  IK ) ) . GT . 0 . 0 . AND •  AHS ( CEL ( K ) ) . LT . 0 . 5 )  FAC-0.0 
IF  (ARS (CEL (K) ) .GF. 0.6. AMD.  AHS (CEL <K) ) .LT.0.9)  FAC-O. 2 
IF  (ARS(CFL(K) ) .GE.O.R.AND.  AHS (CEL (K ) ) »LT.  1 .?)  FAC-0.3 
IF  (ARS(CEL(K)) .GE. 1.2. AND.  AHS(CELIK) ) .LT.1.5)  FAC-0.4 
IF  (ARStCELIK) ) .GE.l.S.AND.  AHS (CEL (K) ) .LT. 1 .9)  FAC-0.5 
IF  (AHS(CEL(K) ) .GE. 1.9. AMD.  AHS (CEL (K) ) .LT.2.2)  FAC-0.6 
IF  (  ARC  (CEL  <K)  )  .GE.2.2.AN0.  A8S  ( CEL  ( K  )  )  .  LT  ,  2  .  *» )  FAC-O. 7 
IF  (AHS(CEL(K) ) .GE.2.S.AND.  A8S (CEL (K 1 ) .LT.2.9)  FAC-O. 8 
IF  (ARS (CEL (K) ) . GE.2.9.AND.  AHS ( CEL ( K ) ) .LT , 3. 2)  FAC-0.9 
IF  (A«S(CEL(K) ) .GE. 3.2. AMD.  AHS (CEL  IK > > .LT . 3. b )  FAC-1.0 


LITS  80S 
LITS  807 
LITS  BOR 
LITS  R 09 
LITS  810 
LTTS  «11 
LITS  b’2 
LTTS  813 
LTTS  814 
LITS  HIS 

l rrs  bis 

LITS  017 
LITS  018 
LTTS  819 
LITS  820 
LITS  821 
LITS  82? 
LITS  823 
LITS  024 
LITS  S?S 


IF 

(ARS (CEL (K) ) 

.  GF  . 

1.2 

.AND. 

AHS ( CEL ( K ) ) 

.LT. 

1.6) 

FAC  = 

.4 

LITS 

8?6 

IF 

(ARS (CEL  (K) ) 

.GE. 

1.6 

.ANO. 

AHS (CEL  <K ) ) 

.LT. 

1.9) 

FAC  = 

.S 

LITS 

827 

IF 

(ARStCELIK) ) 

.  GE  • 

1.9 

.AND. 

AHS (CEL (K) ) 

.LT. 

2.3) 

FAC  = 

.6 

LITS 

828 

IF 

(AHSICEL (K ) ) 

.  GE  • 

2.3 

.AND. 

AHS ( CEL ( K ) ) 

.LT. 

2.6) 

EAC  = 

.7 

LITS 

829 

IF 

( AHS ( CEL ( K ) ) 

.GE. 

2.6 

.AND. 

AHSICEL ( K ) ) 

.LT. 

2.9) 

EAC  = 

•  H 

LITS 

630 

IF 

( ARS ( CEL ( K ) ) 

.GE. 

2.9 

.AND. 

4BS (CEL  (K)  ) 

.LT. 

3.3) 

FAC* 

.9 

LTTS 

831 

IF 

(ARSICElOd  ) 

.GE. 

3.3 

.AND. 

ABS (CEL  (K)  ) 

.LT. 

3.6) 

F  aC  =  1 

.0 

LITS 

83? 

IF 

( AHS ( CEL (K> ) 

.GE. 

3.6 

.ANO. 

ARS ( CEL  <  k ) ) 

.LT. 

4.0) 

F  AC- 1 

.1 

LITS 

833 

IF 

(ARS (CEL  (K)  ) 

.  GE  • 

4.0 

.AND. 

ABS (CEL  (K  )  ) 

.LT. 

4.3) 

FAC-  1 

.2 

LITS 

034 

IF 

( ARS ( CEL  IK) ) 

.GE. 

4.3 

.AND. 

ABS (CEL  (K  )  ) 

.LT. 

4.6) 

F  AC-  1 

.3 

LITS 

835 

GO 

TO  4 

LTTS 

836 

LITS  837 
LITS  838 
LITS  B39 
LITS  040 
LITS  841 
LITS  84? 
LITS  843 
LITS  844 
LITS  84S 
LITS  846 
UTS  B4T 
LITS  848 


1C  t  AH^  ICtL  <6  )  )  .6*  .  .1. ->.  «wu.  AHblLtllM).LI.-U«l  fAL  =  UI  LllS 

IF  (ARSTCEL (K) > .GE.3.H.AN0.  AdS < CEL ( K ) ) .L T . 4 . 2 )  FAC=1.2  LITS 

IF  <ARS(CEL<*> ) .GE.4.2.ANL).  ABS (CEL (K ) ) .LT.4.S)  FAC=1.3  LITS 

IF  (ABS(CELIK) ) .GE.4.S.AM0.  ABS (CEL ( K ) )  .  LT.4.8)  FAC  =  1.4  LITS 

IF  (  ABS (CEL (K) ) .GE. 4. B.ANO.  ABS  (CEL  (K)  )  .LT.5.2)  FAC=1.5  LTTS 

:  FINO  DESIRED  LOCATION  wITh  CURRENT  SPEEDS  AND  DIRECTIONS  FOR  LITS 

:  12  HOUR  PERIOD  LITS 

14  Rfc AD ( 1 ,  103) L AD.LAM.LOO.LOM, <  CUR ( I ) . MO  I R ( I) . I  =  1  *  1 2 ) . I A W  LITS 

103  F0PMAT(I2.I?,1X«I3»14,1<«1?(F?.1%I2»1X) »SX.T2)  I.TTS 

IF  (LOO  .EQ.  999)  00  TO  93  LITS 

IF ( 37. EQ. ILAOJGO  TO  104  LITS 

104  IF  (  ILAH.6F.4H.  AND.  ILAM.LE.S0)  GO  TO  1  OS  I.TTS 

SO  TO  114  LITS 

105  IF(ILOM.(5E.22So.AND.ILOH.LT.3?SO)GO  TO  106  LITs 

00  TO  114  LITS 

106  KLOM=LOM ♦ 126  LTTS 

KLAM=LAM*1  LTTS 

IF (LOn.EO. ILOD) GO  TO  136  LITS 

GO  TO  14  LI  TO 

136  IF (ILOM.GE.LOM.AND. I  LOM. LT . <LOM > GO  TO  137  LTTS 

GO  TO  14  LTTc 

137  IF  (LAO.EO.  ILAD)  GO  TO  139  I.TTS 

GO  TO  14  LTTS 

13R  IF ( ILAH.GE.LAM. AND. ILAm.LT.KLam) 00  TO  2S  LTTS 

GO  TO  14  LTTS 

114  LL0M=L0m*2S"  LTTS 

LL AM=L AM*?  LITS 

IF (LOD.EO. ILOD) GO  TO  116  LTTS 

GO  TO  14  LTTS 

116  IF ( IL0M.GE.L0M.AN0. ILOM.LT.LLOM) GO  TO  117  LITS 

GO  TO  14  LTTS 

117  IF(LAD.EO.ILAO)  GO  TO  119  LTTS 

GO  TO  14  I.ITS 

119  IF(ILAM.GE.LAM.AND.ILAM.LT.LLAM)  00  TO  2S  LITS 

GO  TO  14  LTTS 

:  FIND  DESIRED  LOCATION  with  CURRENT  SPEEDS  AND  DIRECTIONS  FOR  LITS 

:  13  HOUR  PERIOD  IN  LONG  ISLAND  SOUND  LITS 

4  REAO( 1 ,3)LAD»LAM,L0D.L0M, (CUR ( I ) .MOIR ( I ) . 1=1 » 1 3) . I AR  LTTS 

3  F0RMAT(I2,I2,I3.I2.1X,13(F2.1,I2,1X) .12)  LTTS 

IF  (L0r>.E0.999)G0  TO  9*  LITS 

JL0M=L0M*3  LTTS 

JLAM=LAM*?  LTTS 

IF  (Lon  .EO.  ILOD)  GO  TO  6  LITS 

RO  TO  4  LTTS 

6  IF  (ILOM  .GE.  LOM  .aNii.  IlOM  .LT.  JLOM)  GO  TO  7  LITS 

GO  TO  4  LTTS 

7  IF  (LAO  .EO.  ILAD)  GO  TO  9  LTTS 

GO  TO  4  LTTS 

O  IF  (ILAM  .GE.  lam  .AND.  TLAM  .LT.  JLAM)  GO  TO  4025  LITS 

GO  TO  4  LTTS 

4025  WRITF  (61,4026)  L  I  T9 

4026  FORMAT  (1*0.15X.*CUHWENT  SPEEDS  FOR  E AC*  HOUR  OF  T*E  CURRENT*)  LITS 

WRITE(61,40?7)  (CUR(I) .1=1,13)  LITS 

4027  FORMAT ( 1  mo » 1 3 ( F4 , 1 » 2* ) )  LITS 

X JJ1 JJ1 /PERCOL  LTTS 

XJJ1»XJJ1*.S  LITS 

JJ1*XJJ1  LTTS 

IF  (CFL(K)  .LT.  0.)  GO  TO  4040  LTTS 

IF  (JJl  .NF.  o)  JJ1=JJ1*1  LTTS 

IF  (JJI.GT.6)  JJ1=6  LTTS 

IF  (JJ1  .EO.  0)  JJ1=1  LITS 

IF  ( JSLAX (K) .F0.9999)  GO  TO  4099  LTTS 

IF  <JSL«K ( K ♦ 1 ) .EQ.9999)  GO  TO  4099  LITS 

4099  IF  (CUR(JJl)  .EO.  0.0)  GO  TO  100  LITS 

r cuy»ruw  t  jji : lit* 

NOIR*MDIR(JJ1)*10  li." 


(HI  ID  4(1*10 

L  11  S 

VIS 

4040 

J  J?  =  J  J  1*7 

LITS 

9)4 

IF  <  J]?.fiT.l  3)  J  J?=  1  3 

LTTS 

91  7 

IF  (ChR(JJ2)  .til.  0.0)  C.D  To  100 

LITS 

919 

FCUR=CUR ( Jj?) *EAC 

LITS 

919 

NOI«=mOIR( JJ2) *10 

LITS 

920 

40S0 

K  =KK< 

LTTS 

921 

IF  (IDISC.EQ.?)  00  TO  -I/' 

LTTS 

9?? 

WRITE  (M  ,4051)  MO.NriDY.IYwl.MTiME.JSLAK (K) ,FCuR*NDIR* ILAD, 

ILAM, 

LITS 

927 

1  ILOO.  ILOM.FAC 

LITS 

924 

405  1 

FOOVAT  (  1H0  *  3flOlM  ,  I 2  *  1 H/ ♦ I ? .  1 H/ ,  I 4 , 4H  AT  .14.".  WITH  SLACK 

WATER 

ATLITS 

925 

1  THF  PACE  AT  ». 14. 1H, ,/.»  THE  CURRENT  VELOCITY  IS  ».F4. 

2,10 

LTTS 

9?4 

?h  KNOTS  AT  .13,*  DEGREES  aT  POSITION  * , I 2* 1  *  » I ?*2hN, , 1 X 

,13 

. IX.  12 

,?lits 

927 

3  H  w  •  *  factor  IS  *«F3.1) 

LITS 

920 

no  TO  9.3 

LTTS 

929 

PS 

WRITE  (61,26) 

LITS 

930 

26 

FOP HAT  ( 1H0, 15<,*CURREnT  SPEEDS  FOP  Each  HOOP  OF  THE  CURRENT*) 

LTTS 

931 

WRITE (41 , ?  7 )  (CUP ( I ) , 1  =  1 , 12) 

LITS 

97? 

27 

FORMAT (1ho,12(F4.1,?X) ) 

LITS 

933 

XJJ1=XJJ1/PERC0L 

LITS 

934 

x JJl=x JJl ♦ .5 

LTTS 

975 

JJ1=X JJl 

LTTS 

934 

IF  (CFL(K)  .LT.  0.)  GO  TO  40 

LITS 

977 

IF  (J.Jl  .MU.  0)  JJ1=JJ1*1 

LITS 

97« 

I F  (JJl.GT.h)  JJl  =4 

LTTS 

939 

IF  (J.ll  •  F  (2  ,  0)  JJ1=1 

LTTS 

940 

IE  (JSLAK  (  K  !  •  F  Cl .  9994  )  GO  TO  49 

LITS 

941 

IF  (  JSLAK  (K*  1  )  .to. 9444)  Of)  TO  94 

LTTS 

94? 

04 

IF  (0 1°  (  J  J 1  1  .  Fd  .  0.0)  GO  TO  100 

LITS 

947 

FCUP=CUH ( JJ1 ) *F AC 

LITS 

944 

Mi)IR  =  H0T9  (JJ1  )  *10 

LTTS 

945 

G(>  TO  SO 

LTTS 

944 

40 

JJ?=JJ1 ♦ 7 

LTTS 

947 

IF  (JJ2.GT.12)  J  J2=  1  2 

LFTS 

948 

IF  <CHR(JJ?>  . Ed •  0.0)  GO  TO  100 

LITS 

949 

FCUP  =  CUP ( JJ?) *F  AC 

LITS 

950 

N0IP=HPIR(JJ2) *10 

LITS 

951 

SO 

K=KKK 

LTTS 

95? 

IF  (IOTSC.FO.P)  GO  TO  S2 

LITS 

953 

WHITE (M .SI )  “O.NBOY. I YP1  .NT  IMF .JSLAK (K)  . F CUP , NO  I R , I L AD 

,  I LAM, 

LITS 

954 

1 ILOO. ILOM.FAC 

LITS 

955 

si 

FORMAT ( 1H0.3HON  ,  I?,  1H/, 12, In/. I4.4H  aT  ,I4.»,  WITH  SLACK 

WATER 

ATLITS 

956 

1  THE  D ACE  AT  ",  14, 1H, ,/,»  THE  CJRRENT  VELOCITY  IS  »,F4. 

2,10 

LITS 

957 

?h  KNOTS  AT  ,17,*  DEGREES  a  T  POSITION  *, 12. 1 X , I2.2HN. *  1 X 

.13 

.IX. 12 

.2LITS 

988 

3Hw.*  facto0  IS  * , F 3 . 1 ) 

LTTS 

959 

GO  TO  93 

LTTS 

940 

5? 

WPI TE (M .S3)  M0 , NHD Y , I YR 1 .NT  IMF. JSLAK (K)  . FCUR . NO  I R . I L AD 

•  I L  AM  « 

LITS 

941 

1  ILOO, ILOM, FAC 

LITS 

94? 

53 

FORMAT ( 1HO,3HON  ,  12, 1H/, 12, 1H/, 14, 4h  AT  ,I4.»,  WITH  SLACK 

WATER 

ATLITS 

967 

l  THE  golden  gate  AT  ",  14, 1H. ,/,"  THE  CURRENT  VELOCITY  IS  " 

*F 4.2, 

10LITS 

964 

2H  KNOTS  AT  ,13,*  DEGREES  AT  POSITION  «  ,  I  2 ,  1  X  ,  1 2,  ?HN.  ,  1  X 

.13 

,  1  X  »  1 4 

»?L I TS 

965 

3Hw . *  FACTOR  IS  **F3. 1 ) 

LTTS 

966 

GO  TO  93 

LITS 

967 

]00 

WRITE  (Si, 101)  IL  AO,  I  LAM,  ILOO,  ILOM,  NT  IMF. 

LTTS 

968 

101 

FORMAT ( IMA, *TMtHE  IS  NO  CURRENT  AT  POSITION  *,/. 

LTTS 

969 

U3,1X.I2.?hn,,1X,I3,1X,I4,1h-,*  AT  * , I  4 ) 

LITS 

970 

REWIND  l 

LITS 

971 

RETURN 

LITS 

97? 

94 

nnpay  =  0 

LITS 

973 

RETURN 

LTTS 

974 

F.ND 

LITS 

978 

FIMIC 
1i*nFF  <C.  ,LITS) 

S*OFF<m,,C<iOii-Tam.lTTSF.01  .CGAS.DmO?...  .  .  ,  I .  999999 ) 
t*r>EF  <w,  ,CGO'i-TAM.LTTSF.02*CGAG,UP0?.UNIISEO) 

•  • 
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